TITLE OF THE INVENTION 



SINGLE NUCLEOTIDE POLYMORPHISMS 
AND THEIR USE IN GENETIC ANALYSIS 



FIELD OF THE INVENTION 

The present invention is in the field of recombinant DMA 
technology. More specifically, the invention is directed to 
molecules and methods suitable for identifying single nucleotide 
polymorphisms in the genome of an animal, especially a horse or a 
human, and using such sites to analyze identity, ancestry or genetic 
traits. 

CROSS-REFERENCE TO RELATED APPLICATIONS 

l«J^ (\^\ 

This application is^a continuation-in-part of U.S. Patent 
Application Serial No. 08/145,145 (filed November 3, 1993).'V^ /^-Q 
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BACKGROUND OF THE INVENTION 

The capacity to genotype an animal, plant or microbe is of 
fundamental importance to forensic science, medicine and 
epidemiology and public health, and to the breeding and exhibition 
of animals. Such a capacity is needed, for example, to determine 
the identity of the causative agent of an infectious disease, to 
determine whether two individuals are related, or to establish 
whether a particular animal such as a horse is a thoroughbred. 

The analysis of identity and parentage, along with the 
capacity to diagnose disease is also of central concern to human, 
animal and plant genetic studies, particularly forensic or paternity 
evaluations, and in the evaluation of an individual's risk of genetic 
disease. Such goals have been pursued by analyzing variations in 



DNA sequences that distinguish the DNA of one individual from 
another. 

If such a variation alters the lengths of the fragments that 
are generated by restriction endonuclease cleavage, the variations 
are referred to as restriction fragment length polymorphisms 
("RFLPs"). RFLPs have been widely used in human and animal 
genetic analyses (Glassberg, J., UK patent Application 2135774; 
Skolnick, M.H. et al.. Cvtoaen. Cell f^PnPt :^9-'^R.f^7 (1982); Botstein! 
D. et_aL, Ann. J. Hum. Genet. 32:314-331 (1980); Fischer, S.G et al.' 
(PCT Application WO90/13668); Uhlen, M., PCT Application 
WO90/11369)). Where a heritable trait can be linked to a 
particular RFLP, the presence of the RFLP in a target animal can be 
used to predict the likelihood that the animal will also exhibit the 
trait. Statistical methods have been developed to permit the 
multilocus analysis of RFLPs such that complex traits that are 
dependent upon multiple alleles can be mapped (Lander, S. et al.. 
Proc. Natl. Acad. Sci. AJ R A \ s.q-7':t^':t.7ctg7 (igss); Lander, S. et al.! 
Proc. Natl. Acad. Sci. (U.S.A.) 84:2363-2367 (1987); Donis-Keller, H. 
et_aL.n6n^:31 9-337 (1987); Lander, S. eLaL, Genetics 111:1 85- 
199 (1989), all herein incorporated by reference). Such methods 
can be used to develop a genetic map, as well as to develop plants 
or animals having more desirable traits (Donis-Keller, H. et al. . Cell 
5J_:319-337 (1987); Lander, S. et_aL, Genetics 12J_: 1 8 5 - 1 9 9 
(1989)). 

In some cases, the DNA sequence variations are in regions of 
the genome that are characterized by short tandem repeats (STRs) 
that include tandem di- or tri-nucleotide repeated motifs of 
nucleotides. These tandem repeats are also referred to as "variable 
number tandem repeat" ("VNTR") polymorphisms. VNTRs have been 
used in identity and paternity analysis (Weber, J.L., U.S. Patent 
5,075,217; Armour, J.A.L et__al.. FEBS Lett. 307:1 1.^-1 l.q (1992); 
Jones, L eLal.. Eur. J. Haematol .qQ-iAA.iAy (1987); Horn, G.T. 
aL. PCT Application WO91/14003; Jeffreys, A.J., European Patent 
Application 370,719; Jeffreys, A.J., U.S. Patent 5,175,082); 
Jeffreys. A.J. et al.. Amer. J. Hum fipnAt ^t:i-M.0A (1986); Jeffreys. 
A.J. el^, Nature ai£: 76-79 (1985); Gray, I.C. et al.. Proc. R. Anari 



Soc. Lond. 24^:241-253 (1991); Moore, S.S. et_aL, Genomms 10:654- 
660 (1991); Jeffreys, A.J. et al. . Anim. Genet. 1 8 :1-15 (1987); 
Hillel, J. et al. . Anim. Genet. 20:145-155 (1989); Hillel, J. et al. . 
Genet. 124 :783-789 (1990)) and are now being used in a large 
nunfiber of genetic mapping studies. 

A third class of DNA sequence variation results from single 
nucleotide polymorphisms (SNPs) that exist between individuals of 
the same species. Such polymorphisms are far more frequent than 
RFLPs, STRs and VNTRs. In some cases, such polymorphisms 
comprise mutations that are the determinative characteristic in a 
genetic disease. Indeed, such mutations may affect a single 
nucleotide in a protein-encoding gene In a manner sufficient to 
actually cause the disease (i.e. hemophilia, sickle-cell anemia, 
etc.). In many cases, these SNPs are in noncoding regions of a 
genome. Despite the central importance of such polymorphisms in 
modern genetics, no practical method has been developed that 
permits the use of highly parallel analysis of many SNP alleles in 
two or more individuals in genetic analysis. 

The present invention provides such an improved method. 
Indeed, the present invention provides methods and gene sequences 
that permit the genetic analysis of identity and parentage, and the 
diagnosis of disease by discerning the variation of single 
nucleotide polymorphisms. 

SUMMARY OF THE INVENTION 

The present invention is directed to molecules that comprise 
single nucleotide polymorphisms (SNPs) that are present in 
mammalian DNA, and in particular, to equine and human genomic 
DNA polymorphisms. The invention is directed to methods for (i) 
identifying novel single nucleotide polymorphisms (ii) methods for 
the repeated analysis and testing of these SNPs in different 
samples and (iii) methods for exploiting the existence of such sites 
in the genetic analysis of single animals and populations of 
animals. 



The analysis (genotyping) of such sites is useful in 
determining identity, ancestry, predisposition to genetic disease, 
the presence or absence of a desired trait, etc. in detail, the 
invention provides a nucleic acid primer molecule having a 
polynucleotide sequence complementary to an "invariant" 
nucleotide sequence of a genomic DNA segment of a mammal, the 
genomic segment being located immediately 3'-distal to a single 
nucleotide polymorphic site, X, of a single nucleotide polymorphic 
allele of the mammal; and wherein template-dependent extension 
of the nucleic acid primer molecule by a single nucleotide extends 
the primer molecule by a single nucleotide, the single nucleotide 
being complementary to the nucleotide, X, of the single nucleotide 
polymorphic allele. The invention particularly concerns the 
embodiment wherein the mammal is selected from the group 
consisting of humans, non-human primates, dogs, cats, cattle, 
sheep, and horses. 

The invention particularly concerns the embodiments wherein 
the mammal is a horse, and wherein the nucleic acid molecule has a 
nucleotide sequence selected from the group consisting of SEQ ID 
N0:(2n+1) [refer to Table 1], wherein n is an integer selected from 
the group consisting of 0 through 35, or wherein the sequence of 
the immediately 3'-distal segment includes a sequence selected 
from the group consisting of SEQ ID NO:(2n+2), wherein n is an 
integer selected from the group consisting of 0 through 35. 

The invention also provides a nucleic acid molecule having a 
sequence complementary to a sequence selected from the group 
consisting of SEQ ID N0:1 through SEQ ID NO:72. The invention also 
provides a set of at least two of such nucleic acid molecules. 

The invention also provides a set of at least two nucleic acid 
molecules, wherein at least one of the nucleic acid molecules has a 
sequence complementary to a sequence selected from the group 
consisting of SEQ ID N0:1 through SEQ ID NO:72. 

The invention also provides a method for determining the 
extent of genetic similarity between DNA of a target horse and DNA 
of a reference horse, which comprises the steps: 



A) determining, for a single nucleotide polymorphism of the 
target horse, and for a corresponding single nucleotide 
polymorphism of the reference horse, whether the 
polymorphisms contain the same single nucleotide at their 
respective polymorphic sites; and 

B) using the comparison to determine the extent of genetic 
similarity between the target horse and the reference 
horse. 

The invention also concerns the embodiment of such method 
wherein the polymorphic sites are flanked by (1) an immediately 
5'-proximal sequence selected from the group consisting of SEQ ID 
N0:(2n+1), and (2) an immediately 3'-distal sequence selected from 
the group consisting of SEQ ID NO:(2n+2); wherein n is an integer 
selected from the group consisting of 0 through 35. 

The invention particularly concerns the embodiment wherein, 
in step A, the determination Is accomplished by a method having 
the sub-steps: 

(a) incubating a sample of nucleic acid containing the single 
nucleotide polymorphism of the target horse, or the single 
nucleotide polymorphism of the reference horse, in the 
presence of a nucleic acid primer and at least one 
dideoxynucleotide derivative, under conditions sufficient 
to permit a polymerase mediated, template-dependent 
extension of the primer, the extension causing the 
incorporation of a single dideoxynucleotide to the 3'- 
terminus of the primer, the single dideoxynucleotide being 
complementary to the single nucleotide of the 
polymorphic site of the polymorphism; 

(b) permitting the template-dependent extension of the 
primer molecule, and the incorporation of the single 
.dideoxynucleotide; and 

(c) determining the identity of the nucleotide incorporated 
into the polymorphic site, the identified nucleotide being 
complimentary to the nucleotide of the polymorphic site. 

The invention further concerns the embodiment of the above 
methods wherein the template-dependent extension of the primer 



is conducted in the presence of at least two dideoxynucleotide 
triphosphate derivatives selected from the group consisting of 
ddATP, ddTTP, ddCTP and ddGTP, but in the absence of dATP, dTTP, 
dCTP and dGTP. 

The invention particularly concerns the sub-embodiments of 
the above methods wherein the nucleic acid of the sample is 
amplified in vitro prior to the incubation, and/or the primer is 
immobilized to a solid support. 

The invention further concerns the embodiment of the above 
methods wherein a non-invasive swab is used to collect the sample 
of DNA. 

The invention further provides a method for determining the 
probability that a target horse will have a particular trait, which 
comprises the steps: 

A) determining the identity of a single nucleotide present at 
a polymorphic site of an equine single nucleotide 
polymorphism, and being present in more than 51% of a 
set of reference horses; 

B) determining whether a single nucleotide present at a 
polymorphic site of a corresponding single nucleotide 
polymorphism of the target horse has the same identity as 
the single nucleotide present at the polymorphic site of 
the 51% of reference horses exhibiting the trait; 

C) using the determination of step B to establish the 
probability that the target horse will have the particular 
trait. 

The invention further provides a method for creating a 
genetic map of unique sequence equine polymorphisms which 
comprises the steps: 

A) identifying at least one pair of inter-breeding reference 
horses, wherein each of the pairs of horses is characterized 
by having a first and a second reference horse, 
the first reference horse having: 

two alleles (i) and (ii), the alleles each being single 
nucleotide polymorphic alleles having a single 
nucleotide polymorphic site; 



the second reference horse having: 

a corresponding allele (i') to the allele (i) of the 
first reference horse, wherein the allele (i') has a 
single nucleotide polymorphic site, and wherein the 
single nucleotide present at the polymorphic site of 
the allele (i") differs from the single nucleotide 
present at the polymorphic site of the allele (i) of 
the first reference horse, and 

B) identifying in a progeny of at least one of the pairs of inter- 
breeding reference horses the single nucleotide present at a 
single nucleotide polymorphic site of a corresponding allele 
of the alleles (i) and (i'), and the single nucleotide present at 
a single nucleotide polymorphic site of a corresponding 
allele of the alleles (ii) and (ii'); and 

C) determining the extent of genetic linkage between the 
alleles (i) and (ii), to thereby create the genetic map. 

The invention further provides a method for predicting 
whether a target horse will exhibit a predetermined trait which 
comprises the steps: 

A) identifying one or more alleles associated with the trait, 
each allele being a single nucleotide polymorphic allele 
having a single nucleotide polymorphic site; 

B) determining for each of the single nucleotide polymorphic 
alleles, a nucleotide present at the allele's polymorphic 
site in a reference horse exhibiting the trait, to thereby 
define a set of single nucleotides at a set of polymorphic 
sites that are present in a reference horse exhibiting the 
trait; 

C) determining the identity of single nucleotides present at 
corresponding single nucleotide polymorphic alleles of the 
target horse; and 

D) comparing the identity of the single nucleotides present at 
the polymorphic sites of the polymorphisms of the 
reference animal with the single nucleotides present at the 
corresponding single nucleotide polymorphic alleles of the 
target horse. 



The invention further provides a method for identifying a 
ngle nucleotide polymorphic site which comprises: 

A) isolating a fragment of genomic DNA of a reference 
organism; 

B) sequencing the fragment of DNA to thereby determine the 
nucleotide sequence of a segment of the fragment, the 
segment being of a length sufficient to define the 
nucleotide sequence of a pair of oligonucleotide primers 
capable of mediating the specific amplification of the 
fragment; 

C) using the oligonucleotide primers to mediate the specific 
amplification of DNA obtained from a plurality of other 
organisms of the same species as the reference organism; 
and 

D) determining the nucleotide sequences of the amplified DNA 
molecules of step C, and comparing the sequence of the 
amplified molecules with the sequence of the fragment of 
the reference organism to thereby identify a single 
nucleotide polymorphic site. 

The invention also includes a method for interrogating a 
polymorphic region of a human single nucleotide 
polymorphism of a target human, the method comprising: 

A) selecting a known human single nucleotide polymorphism 
for interrogation; 

B) identifying the sequence of at least one oligonucleotide 
that flanks the selected single nucleotide polymorphism; 
the identified sequence being of a length sufficient to 
permit the identification of primers capable of being used 
to effect the specific amplification of the flanking 
oligonucleotide and the polymorphism; 

C) using the primers to effect the amplification of the 
flanking oligonucleotide and the polymorphism of the 
single nucleotide polymorphism of the target human; and 

D) interrogating the single nucleotide polymorphism of the 
amplified polymorphism by genetic bit analysis. 
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BRIEF DESCRIPTION OF THE FIGURES 

Figure 1 illustrates the preferred method for cloning random 
genomic fragments. Genomic DNA^^ size fractionated, and then 
5 introduced into a plasmid vector, in order to obtain random clones. 
PGR primers are designed, and used to sequence the inserted 
genomic sequences. 

Figure 2 illustrates the data generated by preferred method 
for identifying new polymorphic sequences which is cycle 

1 0 sequencing of a random genomic fragment. 

Figure 3 illustrates the RFLP method for screening random 
clones for polymorphic sequences. After the initial optimization of 
PGR conditions (top panel), amplified material is cleaved with 
several restriction enzymes, and the resulting profiles are 

1 5 analyzed (middle panels). A population study is then performed to 
determine allelic frequencies. 

Figure 4 shows a graph of the probability that two 
individuals will have identical genotypes with given panels of 
genetic markers. The number of tests employed is plotted on the 

20 abscissa while the cumulative probability of non-identity is 
plotted on the ordinate. The horizontal line indicates 0.95 
probability of non-identity. Legend: o indicates the extrapolated 
prototype; x indicates 3 alleles (51%, 34%, 15%); triangle indicates 
2 alleles (79%, 21%). 

25 Figure 5 shows a graph of the probability that given panels 

of 20 genetic markers will exclude a random alleged father in a 
paternity suit in which the mother is not in question. The number 
of tests employed is plotted on the abscissa while the cumulative 
probability of exclusion is plotted on the ordinate. The horizontal 

30 line indicates 0.95 probability of exclusion. The legend is as in 
Figure 4. 

Figure 6 uses the SNP identified in clone 177-2 to illustrate 
the organization of the sequences in Table 1. 

Figure 7 illustrates the preferred method for genotyping 
35 SNPs. The seven steps illustrate how GBA can be performed 
starting with a biological sample. 
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Figures 8A and 8B illustrate how horse parentage data 
appears at the microtiter plate level. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

I. The Single Nucleotide Polymorphisms of the Present 
Invention and The Advantages of their Use in Genetic 
Analysis 

A. The Attributes of the Polymorphisms 

The particular gene sequences of interest to the present 
invention comprise "single nucleotide polymorphisms." A 
"polymorphism" is a variation in the DNA sequence of some 
members of a species. The genomes of animals and plants naturally 
undergo spontaneous mutation in the course of their continuing 
evolution (Gusella, J. P., Ann. Rev. Biochem. 55:831-854 (1986)). 
The majority of such mutations create polymorphisms. The 
mutated sequence and the initial sequence co-exist in the species' 
population. In some instances, such co-existence is in stable or 
quasi-stable equilibrium. In other instances, the mutation confers 
a survival or evolutionary advantage to the species, and 
accordingly, it may eventually (i.e. over evolutionary time) be 
incorporated into the DNA of every member of that species. 

A polymorphism is thus said to be "allelic," in that, due to the 
existence of the polymorphism, some members of a species may 
have the unmutated sequence (i.e. the original "allele") whereas 
other members may have a mutated sequence (i.e. the variant or 
mutant "allele"). In the simplest case, only one mutated sequence 
may exist, and the polymorphism is said to be diallelic. Diallelic 
polymorphisms are the most common and the preferred 
polymorphisms of the present invention. The, occurrence of 
alternative mutations can give rise to / ^itil l o l o i o, etc. 
polymorphisms. An allele may be referred to by the nucleotlde(s) 
that comprise the mutation. Thus, for example, in Table 1, clone 
177-2 (SEQ ID N0:1 and SEQ ID N0:2) illustrates the sequence of one 



strand of a diallelic polymorphism in which one allele has a "C" and 
the other allele has a "T" at the polymorphic site. 

The present invention is directed to a particular class of 
allelic polymorphisms, and to their use in genotyping a plant or 
animal. Such allelic polymorphisms are referred to herein as 
"single nucleotide polymorphisms," or "SNPs." "Single nucleotide 
polymorphisms" are defined by the following attributes. A central 
attribute of such a polymorphism is that it contains a polymorphic 
site, "X," most preferably occupied by a single nucleotide, which is 
the site of variation between allelic sequences. A second 
characteristic of an SNP is that its polymorphic site "X" is 
preferably preceded by and followed by "invariant" sequences of the 
allele. The polymorphic site of the SNP is thus said to lie 
"immediately" 3' to a "S'-proximal" invariant sequence, and 
"immediately" 5* to a "S'-distal" invariant sequence. Such 
sequences flank the polymorphic site. 

As used herein, a sequence is said to be an "invariant" 
sequence of an allele if the sequence does not vary in the 
population of the species, and if mapped, would map to a 
"corresponding" sequence of the same allele in the genome of every 
member of the species population. Two sequences are said to be 
"corresponding" sequences if they are analogs of one another 
obtained from different sources. The gene sequences that encode 
hemoglobin in two humans illustrate "corresponding" allelic 
sequences. The definition of "corresponding alleles" provided 
herein is intended to clarify, but not to alter, the meaning of that 
term as understood by those of ordinary skill in the art. Each row 
of Table 1 shows the identity of the nucleotide of the polymorphic 
site of "corresponding" equine alleles, as well as the invariant S'- 
proximal ^^^^I'p^i^^j^jeqLiences that are also attributes of that 
SNP. ^^yo f rfe jlondionl ^ alleles" are illustrated in Table 5 with 
regard to human alleles. Each row of Table 5 shows the identity of 
the nucleotide of the polymorphic site of "corresponding" human 
alleles, as well as the invariant 5'-proximai and 3'-distal 
sequences that are also attributes of that SNP. 
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Since genomic DNA is double-stranded, each SNP can be 
defined in terms of eitlier strand. Thus, for every SNP, one strand 
will contain an immediately 5'-proximal invariant sequence and the 
other will contain an immediately 3'-distal invariant sequence. In 
5 the preferred embodiment, wherein a SNP's polymorphic site, "X," is 
a single nucleotide, each strand of the double-stranded DNA of the 
SNP will contain both an immediately 5'-proximal invariant 
sequence and an immediately 3'-distai invariant sequence. 

Although the preferred SNPs of the present invention involve 

1 0 a substitution of one nucleotide for another at the SNP's 
polymorphic site, SNPs can also be more complex, and may 
comprise a deletion of a nucleotide from, or an insertion of a 
nucleotide into, one of two corresponding sequences. For example, 
a particular gene sequence may contain an A in a particular 

1 5 polymorphic site in some animals, whereas in other animals a 
single or multiple base deletion might be present at that site. 
Although the preferred SNPs of the present invention have both an 
invariant proximal sequence and invariant distal sequence, SNPs 
may have only an invariant proximal or only an invariant distal 

20 sequence. 

Nucleic acid molecules having -t^^ a sequence complementary 
to that of an immediately 3'-distal invariant sequence of a SNP 
can, if extended in a "template-dependent" manner, form an 
extension product that would contain the SNP's polymorphic site. 

2 5 AST preferred example of such a nucleic acid molecule is a nucleic 

A 

acid molecule whose sequence is the same as that of a 5'-proximal 
invariant sequence of the SNP. "Template-dependent" extension 
refers to the capacity of a polymerase to mediate the extension of 
a primer such that the extended sequence is complementary to the 

3 0 sequence of a nucleic acid template. A "primer" is a single- 

stranded oligonucleotide or a single-stranded polynucleotide that 
is capable of being extended by the covalent addition of a 
nucleotide in a "template-dependent" extension reaction. In order 
to possess such a capability, the primer must have a 3'-hydroxyl 
35 terminus, and be hybridized to a second nucleic acid molecule (i.e. 
the "template"). A primer is typically 11 bases or longer; most 
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preferably, a primer is 20 bases, however, primers of shorter or 
greater length may suffice. A "polymerase" is an enzyme that is 
capable of incorporating nucleoside triphosphates to extend a 3'- 
hydroxyl group of a nucleic acid molecule, if that molecule has 
hybridized to a suitable template nucleic acid molecule. 
Polymerase enzymes are discussed in Watson, J.D., In: Moleniiar 
Biology of the Gene. 3rd Ed., W.A. Benjamin, Inc., MenIo Park, CA 
(1977), which reference is Incorporated herein by reference, and 
similar texts. Other polymerases such as the large proteolytic 
fragment of the DNA polymerase I of the bacterium E. coli. 
commonly known as "Klenow" polymerase, E. coli DNA polymerase I, 
and bacteriophage T7 DNA polymerase, may also be used to perform 
the method described herein. Nucleic acids having the same 
sequence as that of the immediately 3' distal invariant sequence of 
a SNP can be ligated in a template dependent fashion to a primer 
that has the same sequence as that of the immediately 5' proximal 
sequence that has been extended by one nucleotide in a template 
dependent fashion. 

B. The Advantages of Using SNPs in Genetic 
Analysis 

The single nucleotide polymorphic sites of the present 
invention can be used to analyze the DNA of any plant or animal. 
Such sites are particularly suitable for analyzing the genome of 
mammals, including humans, non-human primates, domestic 
animals (such as dogs, cats, etc.), farm animals (such as cattle, 
sheep, etc.) and other economically important animals, in 
particular, horses. They may, however be used with regard to other 
types of animals, particularly birds (such as chickens, turkeys, 
etc.) SNPs have several salient advantages over RFLPs, STRs and 
VNTRs. 

First, SNPs occur at greater frequency (approximately 10- 
100 fold greater), and with greater uniformity than RFLPs and 
VNTRs. The greater frequency of SNPs means that they can be more 
readily identified than the other classes of polymorphisms. The 
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greater uniformity of their distribution permits the identification 
of SNPs "nearer" to a particular trait of interest. The combined 
effect of these two attributes makes SNPs extremely valuable. For 
example, if a particular trait (e.g. predisposition to cancer) 
reflects a mutation at a particular locus, then any polymorphism 
that is linked to the particular locus can be used to predict the 
probability that an individual will be exhibiting that trait. 

The value of such a prediction is determined in part by the 
distance between the polymorphism and the locus. Thus, if the 
locus is located far from any repeated tandem nucleotide sequence 
motifs, VNTR analysis will be of very limited value. Similarly, If 
the locus is far from any detectable RFLP, an RFLP analysis would 
not be accurate. However, since the SNPs of the present invention 
are present approximately once every 300 bases in the mammalian 
genome, and exhibit uniformity of distribution, a SNP can, 
statistically, be found within 150 bases of any particular genetic 
lesion or mutation. Indeed, the particular mutation may itself be 
an SNP. Thus, where such locus has been sequenced, the variation 
in that locus' nucleotide is determinative of the trait in question. 

Second, SNPs are more stable than other classes of 
polymorphisms. Their spontaneous mutation rate is approximately 
10-9, approximately 1,000 times less frequent than VNTRs. 
Significantly, VNTR-type polymorphisms are characterized by high 
mutation rates. 

Third, SNPs have the further advantage that their allelic 
frequency can be inferred from the study of relatively few 
representative samples. These attributes of SNPs permit a much 
higher degree of genetic resolution of identity, paternity exclusion, 
and analysis of an animal's predisposition for a particular genetic 
trait than is possible with either RFLP or VNTR polymorphisms. 

Fourth, SNPs reflect the highest possible definition of 
genetic information - nucleotide position and base identity. 
Despite providing such a high degree of definition, SNPs can be 
detected more readily than either RFLPs or VNTRs, and with greater 
flexibility. Indeed, because DNA is double-stranded, the 
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complimentary strand of the allele can be analyzed to confirm the 
presence and identity of any SNP. 

The flexibility with which an identified SNP can be 
characterized is a salient feature of SNPs. VNTR-type 
polymorphisms, for example, are most easily detected through size 
fractionation methods that can discern a variation in the number of 
the repeats. RFLPs are most easily detected by size fractionation 
methods following restriction digestion. 

In contrast, SNPs can be characterized using any of a variety 
of methods. Such methods include the direct or indirect sequencing 
of the site, the use of restriction enzymes where the respective 
alleles of the site create or destroy a restriction site, the use of 
allele-specific hybridization probes, the use of antibodies that are 
specific for the proteins encoded by the different alleles of the 
polymorphism, or by other biochemical interpretation. 

The "Genetic Bit Analysis {"GBA") method disclosed by Goelet, 
P- et_ai. (WO 92/15712^-d:K3 roin incorporated by rGfcrcno e), and 
discussed below, is a preferred method for detecting the single 
nucleotide polymorphisms of the present invention. GBA is a 
method of polymorphic site interrogation in which the nucleotide 
sequence information surrounding the site of variation in a target 
DNA sequence is used to design an oligonucleotide primer that is 
complementary to the region immediately adjacent to, but not 
including, the variable nucleotide in the target DNA. The target 
DNA template is selected from the biological sample and hybridized 
to the interrogating primer. This primer is extended by a single 
labeled didebxynucleotide using DNA polymerase in the presence of 
two, and preferably all four chain terminating nucleoside 
triphosphate precursors. Cohen, D. et al. (PCT Application 
WO91/02087) describes a related method of genotyping. 

Recently, several primer-guided nucleotide incorporation 
procedures for assaying polymorphic sites in DNA have been 
described (Komher, J. S. et al.. Nucl. Acids. Res. 17:7779-7784 
(1989); Sokolov, B. P., Nucl. Acids Rfi.<; lfi-.?fi7i (1990); Syvanen. A.- 
C, eLal., Genomics 8:684 - 692 (1990); Kuppuswamy, M.N. et al .. 
Proc. Natl. Acad. Sci. Hi .'^ A ) Sft-1 14':>-1 i^t (1991); Prezant, T.R. et 
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al.. Hum. Mutat. 1:159-164 (1992); Ugozzoli, L et al.. GATA 9:107- 
112 (1992); Nyren, P. et al.. Anal. Biochem. 208 :171-175 (1993)). 
These methods differ from GBA in that they all rely on the 
incorporation of labeled deoxynucleotides to discriminate between 
bases at a polymorphic site. In such a format, since the signal is 
proportional to the number of deoxynucleotides incorporated, 
polymorphisms that occur in runs of the same nucleotide can result 
in signals that are proportional to the length of the run (Syvanen, 
A.-C, et al .. Amer. J. Hum. Genet. 52:46-59 (1993)). Such a range of 
locus-specific signals could be more complex to interpret, 
especially for heterozygotes, compared to the simple, ternary (2:0, 
1:1, or 0:2) class of signals produced by the GBA method. In 
addition, for some loci, incorporation of an incorrect 
deoxynucleotide can occur even in the presence of the correct 
dideoxynucleotide (Komher, J. S. et al .. Nucl. Acids. Res. 17:7779- 
7784 (1989)). Such deoxynucleotide misincorporation events may 
be due to the Km of the DNA polymerase for the mispaired deoxy- 
substrate being comparable, in some sequence contexts, to the 
relatively poor Km of even a correctly base paired dideoxy- 
substrate (Kornberg, A., eLai., In: DNA Replication, 2nd Edition, W.H. 
Freeman and Co., (1992); New York; Tabor, S. et al .. Proc. Natl. Acad. 
Sci. m.S.A.t 86:4076-4080 (1989)). This effect would contribute to 
the background noise in the polymorphic site interrogation. 

II. Methods for Discovering Novel Polymorphic Sites 

A preferred method for discovering polymorphic sites 
involves comparative sequencing of genomic DNA fragments from a 
number of haploid genomes. In the preferred embodiment, 
illustrated in Figure 1, such sequencing is performed by preparing a 
random genomic library that contains 0.5-3 kb fragments of DNA 
derived from one member of a species. Sequences of these 
recombinants are then used to facilitate PGR sequencing of a 
number of randomly selected individuals of that species at the 
same genomic loci. 
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From such genomic libraries (typically of approximately 
50,000 clones), several hundred (200-500) individual clones are 
purified, and the sequences of the termini of their inserts are 
determined. Only a small amount of terminal sequence data (100- 
200 bases) need be obtained to permit PGR amplification of the 
cloned region. The purpose of the sequencing is to obtain enough 
sequence information to permit the synthesis of primers suitable 
for mediating the amplification of the equivalent fragments from 
genomic DNA samples of other members of the species. Preferably, 
such sequence determinations are performed using cycle 
sequencing methodology. 

The primers are used to amplify DNA from a panel of 
randomly selected members of the target species. The number of 
members in the panel determines the lowest frequency of the 
polymorphisms that are to be isolated. Thus, if six members are 
evaluated, a polymorphism that exists at a frequency of, for 
example, 0.01 might not be identified. In an illustrative, but 
oversimplified, mathematical treatment, a sampling of six 
members would be expected to identify only those polymorphisms 
that occur at a frequency of greater than about .08 (i.e. 1.0 total 
frequency divided by 6 members divided by 2 alleles per genome). 
Thus, if one desires the identification of less frequent 
polymorphisms, a greater number of panel members must be 
evaluated. 

Cycle sequence analysis (Mullis, K. et al. . Cold Spring Harbor 
Svmp. Quant. Biol. 51 :2fi3-27.'^ (1986); Eriich H. et al.. European 
Patent Appln. 50,424; European Patent Appln. 84,796, European 
Patent Application 258,017, European Patent Appln. 237,362; 
Mullis, K., European Patent Appln. 201,184; Mullis K. et al.. U.S. 
Patent No. 4,683,202; Eriich, H., U.S. Patent No. 4,582,788; and 
Saiki, R. et al.. U.S. Patent No. 4,683,194)) is facilitated through 
the use of automated DNA sequencing instruments and software 
(Applied Biosystems, Inc.). Differences between sequences of 
different animals can thereby be identified and confirmed by 
inspecting the relevant portion of the chromatograms on the 
computer screen. Differences are interpreted to reflect a DNA 
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polymorphism only if the data was available for both strands, and 
present in more than one haploid example among the population of 
animals tested. Figure 2 illustrates the preferred method for 
identifying new polymorphic sequences which is cycle sequencing 
of a random genomic fragment. The PGR fragments from five 
unrelated horses were electroeluted from acrylamide gels and 
sequenced using repetitive cycles of thermostable Taq DNA 
polymerase in the presence of a mixture of dNTPs and fluorescent 
ddNTPs. The products were then separated and analyzed using an 
automated DNA sequencing instrument of Applied Biosystems, Inc. 
The data was analyzed using ABI software. Differences between 
sequences of different animals were identified by the software and 
confirmed by inspecting the relevant portion of the chromatograms 
on the computer screen. Differences are presented as "DNA 
Polymorphisms" only if the data is available for both strands and 
present in more than one haploid example among the five horses 
tested. The top panel shows an "A" homozygote, the middle panel an 
"AT" heterozygote and the bottom panel a "T" homozygote. 

Despite the randomized nature of such a search for 
polymorphisms, such sequencing and comparison of random DNA 
clones is readily able to identify suitable polymorphisms. Indeed, 
with respect to the horse, approximately 1/400 nucleotides 
sequenced by these methods would be discovered as the 
polymorphic site of an SNP. 

The discovery of polymorphic sites can alternatively be 
conducted using the strategy outlined in Figure 3. In this 
embodiment, the DNA sequence polymorphisms are identified by 
comparing the restriction endonuclease cleavage profiles generated 
by a panel of several restriction enzymes on products of the PGR 
reaction from the genomic templates of unrelated members. Most 
preferably, each of the restriction endonucleases used will have 
four base recognition sequences, and will therefore allow a 
desirable number of cuts in the amplified products. 

The restriction digestion patterns obtained from the genomic 
DNAs are preferably compared directly to the patterns obtained 
from PGR products generated using the corresponding plasmid 
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templates. Such a comparison provides an internal control which 
indicates that the amplified sequences from the genomic and 
plasmid DMAs derive from equivalent ici. This control also allows 
identification of primers that fr .aitously amplify repeated 
sequences, or multicopy loci, since these will generate many more 
fragments from the genomic DNA templates than from the plasmid 
templates. 

Ml. Methods for Genotyping the Single Nucleotide 
Polymorphisms of the Present Invention 

Any of a variety of methods can be used to identify the 
polymorphic site, "X," of a single nucleotide polymorphism of the 
present invention. The preferred method of such identification 
involves directly ascertaining the sequence of the polymorphic site 
for each polymorphism being analyzed. This approach is thus 
markedly different from the RFLP method which analyzes patterns 
of bands rather than the specific sequence of a polymorphism. 

A. Sampling Methods 

Nucleic acid specimens may be obtained from an individual of 
the species that is to be analyzed using either "invasive" or "non- 
invasive" sampling means. A sampling means is said to be 
"Invasive" if it involves the collection of nucleic acids from within 
the skin or organs of an animal (including, especially, a murine, a 
human, an ovine, an equine, a bovine, a porcine, a canine, or a feline 
animal). Examples of invasive methods include blood collection, 
semen collection, needle biopsy, pleural aspiration, etc. Examples 
of such methods are discussed by Kim, C.H. et al . (J. Virol. 6fi:.'^R7Q- 
3882 (1992)); Biswas, B. et_ai. (Annals NY Acad. Sci. .qfln-t^flP-t^fig 
(1990)); Biswas, B. et al . ( J. Clin. Minrnhini 9Q-PP9R-99':t':t (1991)). 

In contrast, a "non-invasive" sampling means is one in which 
the nucleic acid molecules are recovered from an internal or 
external surface of the animal. Examples of such "non-invasive" 
sampling means include "swabbing," collection of tears, saliva. 
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urine, fecal material, sweat or perspiration, etc. As used herein-, 
"swabbing" denotes contacting an applicator/collector ("swab") 
containing or comprising an adsorbent material to a surface in a 
manner sufficient to collect surface debris and/or dead or sloughed 
off cells or cellular debris. Such collection may be accomplished 
by swabbing nasal, oral, rectal, vaginal or aural orifices, by 
contacting the skin or tear ducts, by collecting hair follicles, etc. 

Nasal swabs have been used to obtain clinical specimens for 
PGR amplification (Olive, D.M. et al. . J. Gen. Virol. 71 :21 41 -P1 d7 

(1990) ; Wheeler, J.G. et al.. Amer. J. Vet. Rp.q f;p-i7QQ-iflnQ 

(1991) ). The use of hair follicles to identify VNTR polymorphisms 
for paternity testing in horses has been described by Ellegren, H. et 
ai. (Animal Genetics 23:133-142 (1992). The reference states that 
a standardized testing system based on PCR-analyzed 
microsatellite polymorphisms are likely to be an alternative to 
blood typing for paternity testing. 

A preferred swab for the collection of DNA will comprise a 
solid support, at least a portion of which is designed to adsorb 
DNA. The portion designed to adsorb DNA may be of a compressible 
texture, such as a "foam rubber," or the like. Alternatively, it may 
be an adsorptive fibrous composition, such as cotton, polyester, 
nylon, or the like. In yet another embodiment, the portion designed 
to adsorb DNA may be an abrasive material, such as a bristle or 
brush, or having a rough surface. The portion of the swab that is 
designed to adsorb DNA may be a combination of the above textures 
and compositions (such as a compressible brush, etc.). The swab 
will, preferably, be specially formed in a substantially rod-like, 
arrow-like or mushroom-like shape, such that it will have a 
segment that can be held by the collecting individual, and a tip or 
end portion which can be placed into contact with the surface that 
contains the sample DNA that is to be collected. In one 
embodiment, the swab will be provided with a storage chamber, 
such as a plastic or glass tube or cylinder, which may have one 
open end, such as a test-tube. Alternatively, the tube may have 
two open ends, such that after swabbing, the collector can pull on 
one end of the swab so as to cause the other end of the swab to be 
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withdrawn into the tube. In yet another embodiment, the tube may 
have two open ends, such that after swabbing, the tube can be 
converted into a column to assist in the further processing of the 
collected DNA. In one embodiment, the end or ends of the storage 
chamber are self-sealing after swabbing has been accomplished. 

The swab or the storage chamber may contain antimicrobial 
agents at concentrations sufficient to prevent the proliferation of 
microbes (bacteria, yeast, molds, etc.) during subsequent storage 
or handling. 

In one embodiment, the swab or storage chamber will contain 
an chromogenic reagent which reacts to the presence of DNA to 
yield a detectable signal that can be identified at the time of 
sample collection. Most preferably, such a reagent will comprise a 
minimum concentration "open-end point" assay for DNA. Such an 
assay is capable of detecting concentrations of nucleic acids that 
range from the minimum detection level of the assay to the 
maximum assay saturation level of the assay. This saturation level 
is adjustable, and can be increased by decreasing the time of 
reaction. Preferred chromogenic reagents include anti-DNA 
antibodies that are conjugated to enzymes, diaminopimelic acid, 
etc. 

B. Amplification-Based Analysis 

The detection of polymorphic sites in a sample of DNA may be 
facilitated through the use of DNA amplification methods. Such 
methods specifically increase the concentration of sequences that 
span the polymorphic site, or include that site and sequences 
located either distal or proximal to it. Such amplified molecules 
can be readily detected by gel electrophoresis or other means. 

The most preferred method of achieving such amplification 
employs PGR, using primer pairs that are capable of hybridizing to 
the proximal sequences that define a polymorphism in its double- 
stranded form. 

In lieu of PGR, alternative methods, such as the "Ligase Ghain 
Reaction" ("LCR") may be used (Barany, F., Proc. Natl. Acad. Sci. 
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(LL^ 88:189-193 (1991). LCR uses two pairs of oligonucleotide 
probes to exponentially amplify a specific target. The sequences of 
each pair of oligonucleotides^^* selected to permit the pair to 
hybridize to abutting sequences of the same strand of the target 
Such hybridization forms a substrate for a template-dependent 
ligase. As with PGR, the resulting products thus serve as a 
template in subsequent cycles and an exponential amplification of 
the desired sequence is obtained. 

In accordance with the present invention, LCR can be 
performed with oligonucleotides having the proximal and distal 
sequences of the same strand of a polymorphic site. In one 
embodiment, either oligonucleotide will be designed to include the 
actual polymorphic site of the polymorphism. In such an 
embodiment, the reaction conditions are selected such that the 
oligonucleotides can be ligated together only if the target molecule 
either contains or lacks the specific nucleotide that is 
complementary to the polymorphic site present on the 
oligonucleotide. 

In an alternative embodiment, the oligonucleotides will not 
include the polymorphic site, such that when they hybridize to the 
target molecule, a "gap" is created (see, Segev. D., PCT Application 
WO 90/01069). This gap is then "filled" with complementary dNTPs 
(as mediated by DNA polymerase), or by an additional pair of 
oligonucleotides. Thus, at the end of each cycle, each single strand 
has a complement capable of serving as a target during the next 
cycle and exponential amplification of the desired sequence is 
obtained. 

The "Oligonucleotide Ligation Assay" ("OLA") (Landegren, U. et 
aL, Science 241:1077-1080 (1 988)) . shares certain similarities 
with LCR and may also be adapted for use in polymorphic analysis. 
The OLA protocol uses two oligonucleotides which are designed to 
be capable of hybridizing to abutting sequences of a single strand 
of a target. OLA, like LCR, is particularly suited for the detection 
of point mutations. Unlike LCR, however, OLA results in "linear" 
rather than exponential amplification of the target sequence. 
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Nickerson, D.A. eLaL have described a nucleic acid detection 
assay that combines attributes of PGR and OLA (Nickerson, D.A. ei 
aL. Proc. Natl. Acad, 8ci (I I .S A } 87:8923-8927 (1990).' In "this 
method, PGR is used to achieve the exponential amplification of 
target DNA, which is then detected using OLA. In addition to 
requiring multiple, and separate, processing steps, one problem 
associated with such combinations is that they inherit all of the 
problems associated with PGR and OLA. 

Schemes based on ligation of two (or more) oligonucleotides 
in the presence of nucleic acid having the sequence of the resulting 
"di-oligonucleotide", thereby amplifying the di-oligonucleotide, are 
also known (Wu, D.Y. M_ai., Genomics 4:560 (1989)), and may be 
readily adapted to the purposes of the present invention. 

Other known nucleic acid amplification procedures, such as 
transcription-based amplification systems (Maiek, L.T. et al .. U.S. 
Patent 5,130,238; Davey. G. et_al.. European Patent Application 
329,822; Schuster et_ai., U.S. Patent 5,169,766; Miller. H.I. et a! .. 
PGT appln. WO 89/06700; Kwoh, D. et al .. Proc. Natl. Arpd <^h' 
(U.S.A.) M:1173 (1989); Gingeras, T.R. et_ai., PGT application WO 
88/10315)), or isothermal amplification methods (Walker, G.T. el 
aj.. Proc. Natl. Acad. Sci (USA) 89:392-396 (1992)) may also be 
used. 

C. Preparation of Single-Stranded DNA 

The direct analysis of the sequence of an SNP of the present 
invention can be accomplished using either the "dideoxy-mediated 
chain termination method," also known as the "Sanger Method" 
(Sanger, P., ejai.. J. Molec. Biol £4:441 (1975)) or the "chemical 
degradation method," "also known as the "Maxam-Gilbert method" 
(Maxam, A.M., et al.. Proc. Natl. Anari. Sci. ([l R A) jA-'^f^n (1977), 
both references herein incorporated by reference). Methods for 
sequencing DNA using either the dideoxy-mediated method or the 
Maxam-Giibert method are widely known to those of ordinary skill 
in the art. Such methods are, for example, disclosed in Sambrook, 
^i-^- Molecular Gionina. a la boratory Manual. 2nd FHitinn , Gold 




Spring Harbor Press. Cold Spring Harbor, New Yorl< (1989), and in 
Zyskind, J.W., et a!.. Recombinant DNA Laboratory Manual , Academic 
Press, inc. . New Yorl< (1988), both lierein incorporated by reference. 

Where a nucleic acid sample contains double-stranded DNA 
(or RNA), or where a double-stranded nucleic acid amplification 
protocol (such as PGR) has been employed, it is generally desirable 
to conduct such sequence analysis after treating the double- 
stranded molecules so as to obtain a preparation that is enriched 
for, and preferably predominantly, only one of the two strands. 

The simplest method for generating single-stranded DNA 
molecules from double-stranded DNA is denaturation using heat or 
alkalai treatment. 

Single-stranded DNA molecules may also be produced using 
the single-stranded DNA bacteriophage Ml 3 (Messing, J. et al .. Meth. 
Enzvmol. 101 :2Q (1983); see also, Sambrook, J., et al. (In: 
Molecular Cloning: A Laboratory Manual . Cold Spring Harbor 
Laboratory Press, Cold. Spring Harbor, NY (1989)). 

Several alternative methods can be used to generate single- 
stranded DNA molecules. Gyllensten, U. et al. . ( Proc. Natl. Acad. 
Sci. m.S.A.^ 85:7652-765fi (1988) and Mihovilovic, M. et al.. 
( BioTechniques 7(1^ 14 (1989)) describe a method, termed 
"asymmetric PCR," in which the standard "PCR" method is conducted 
using primers that are present in different molar concentrations. 
Higuchi, R.G. et al . ( Nucleic Acids Re.'s. 17:5865 (1985)) exemplifies 
an additional method for generating single-stranded amplification 
products. The method entails phosphorylating the 5'-terminus of 
one strand of a double-stranded amplification product, and then 
permitting a 5' -> 3' exonuclease (such as exonuciease) to 
preferentially degrade the phosphorylated strand. 

Other methods have also exploited the nuclease resistant 
properties of phosphorothioate derivatives in order to generate 
single-stranded DNA molecules (Benkovic et al .. U.S. Patent No. 
4,521,509; June 4, 1985); Sayers, J.R. et al . ( Nucl. Acids Res. 
16:791-802 (1988); Eckstein, F. et al .. Biochemistry 15:1685-1691 
(1976); Ott, J. et al.. Biochemistry 26:fi237-8P41 (1987)). 




A discussion of the relative advantages and disadvantages of 
such methods of producing singie-stranded molecules is provided 
by Nikiforov, T. (U?S. paluii t -afipiicalLon— ^o rial no. Q 8 / 00 5 ,061t 
44e£eip-4t:te orporQtod by rof o fonco) . > . 

Most preferably, such sfngie-siranded molecules will be 
produced using the methods described by Nikiforov, T. (U.S. patent 
application serial no. 08/005,061, herein inoorporatod — by- 
-f^iatance^ In brief, these methods employ nuclease resistant 
nucleotides derivatives, and incorporates such derivatives, by 
chemical synthesis or enzymatic means, into primer molecules, or 
their extension products, in place of naturally occurring 
nucleotides. 

Suitable nucleotide derivatives include derivatives in which 
one or two of the non-bridging oxygens of the phosphate moiety of 
a nucleotide has been replaced with a sulfur-containing group 
(especially a phosphorothioate), an alkyi group (especially a methyl 
or ethyl alkyI group), a nitrogen-containing group (especially an 
amine), and/or a selenium-containing group, etc. 

Phosphorothioate deoxyribonucieotide or ribonucleotide 
derivatives (e.g. a nucleoside 5'-0-1-thiotriphosphate) are the 
most preferred nucleotide^; d;arivativ Any of a variety of 
chemical methods may be used to produce such phosphorothioate 
derivatives (see, for example, Zon, G. et al .. Anti-Canc. Drug Pes. 
g.:539-568 (1991); Kim, S.G. et al .. Biochem. BioDhvs. Res. Commun. 
179:1614-1619 (1991); Vu, H. et a! .. Tetrahedron Lett. 32:3005- 
3008 (1991); Taylor, J.W. et al .. Nucl. Acids Res. 1 3:8749-8764 
(1985); Eckstein, F. et al .. Biochemistrv 15:1685-1691 (1976); Ott, 
J. et a!.. Biochemistrv 26 :8237-8241 (1987); Ludwig, J. et al .. J. 
Org. Chem. ^:631-635 (1989), all herein incorporated by 
reference). Phosphorothioate nucleotide derivatives can also be 
obtained commercially from Amersham or Pharmacia. 

Importantly, the selected nucleotide derivative must be 
suitable for in vitro primer-mediated extension and provide 
nuclease resistance to the region of the nucleic acid molecule in 
which it is incorporated. In the most preferred embodiment, it 
must confer resistance to exonucleases that attack double- 
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stranded DNA from the 5'-end (5'->3' exonucleases). Examples of 
such exonucleases include bacteriophage T7 gene 6 exonuclease 
("T7 exonuclease) and the bacteriophage lambda exonuclease ("A, 
exonuclease"). Both 11 exonuclease and X exonuclease are inhibited 
to a significant degree by the presence of phosphorothioate bonds 
so as to allow the selective degradation of one of the strands. 
However, any double-strand specific, 5'->3' exonuclease can be used 
for this process, provided that its activity is affected by the 
presence of the bonds of the nuclease resistant nucleotide 
derivatives. The preferred enzyme when using phosphorothioate 
derivatives is the T7 gene 6 exonuclease, which shows maximal 
enzymatic activity in the same buffer used for many DNA dependent 
polymerase buffers including Taq polymerase. The 5'-^ 3' 
exonuclease resistant properties of phosphorothioate derivative- 
containing DNA molecules are discussed, for example, in Kunkel, 
T.A. (In: Nucleic Acids and Molecular Biolocv . Vol. 2, 124-135 
(Eckstein, F. et al .. eds.), Springer-Verlag, Berlin, (1988)). The 
3'_^5' exonuclease resistant properties of phosphorothioate 
nucleotide containing nucleic acid molecules are disclosed in 
Putney, S.D., et al. ( Proc. Natl. Anad. Sci. (U.S.A.) 78:7350-7354 
(1981)) and Gupta, A.P., et al. ( Nucl. Acids. Res. . 1 2:5897-591 1 
(1984)). 

In addition to being resistant to such exonucleases, nucleic 
acid molecules that contain phosphorothioate derivatives at 
restriction endonuclease cleavage recognition sites are resistant 
to such cleavage. Taylor, J.W., et al. r Nucl. Acids Res. . 13:8749- 
8764 (1985)) discusses the endonuclease resistant properties of 
phosphorothioate nucleotide containing nucleic acid molecules. 

The nuclease resistance of phosphorothioate bonds has been 
utilized in a DNA amplification protocol (Walker, T.G. et al. (Prpc, 
Natl. Acad. Sci. (U.S.A.I 89:392-396 (1992)). In the Walker §i_ai. 
method, phosphorothioate nucleotide derivatives are installed 
within a restriction endonuclease recognition site in one strand of 
a double-stranded DNA molecule. The presence of the 
phosphorothioate nucleotide derivatives protects that strand from 
cleavage, and thus results in the nicking of the unprotected strand 
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by the restriction endonuclease. Amplification is accomplislned by 
cycling the nicking and polymerization of the strands. 

Similarly, this resistance to nuclease attack has been used as 
the basis for a modified "Sanger" sequencing method (Labeit, S. £t 
ai. (DNA 5:173-177 (1986)). In the Labeit et al . method, 35s. 
labeled phosphorothioate nucleotide derivatives were employed in 
lieu of the dideoxy nucleotides of the "Sanger" method. 

In the most preferred embodiment, the phosphorothioate 
derivative is included in the primer. The nucleotide derivative may 
be incorporated into any position of the primer, but will preferably 
be incorporated at the 5'-terminus of the primer, most preferably 
adjacent to one another. Preferably, the primer molecules will be 
approximately 25 nucleotides in length, and contain from about 4% 
to about 100%, and more preferably from about 4% to about 40%, 
and most preferably about 16%, phosphorothioate residues (as 
compared to total residues). The nucleotides may be incorporated 
into any position of the primer, and may be adjacent to one another, 
or interspersed across all or part of the primer. 

In one embodiment, the present invention can be used in 
concert with an amplification protocol, for example, PGR. In this 
embodiment, it is preferred to limit the number of 
phosphorothioate bonds of the primers to about 10 (or 
approximately half of the length of the primers), so that the 
primers can be used in a PGR reaction without any changes to the 
PGR protocol that has been established for non-modified primers. 
When the primers contain more phosphorothioate bonds, the PGR 
conditions may require adjustment, especially of the annealing 
temperature, in order to optimize the reaction. 

The incorporation of such nucleotide derivatives into DNA or 
RNA can be accomplished enzymatically, using a DNA polymerase 
(Vosberg, H.P. et al .. Biochemistrv 16 : 3633-3640 (1977); Burgers, 
P.M.J, et al .. J. Biol. Ghem. 254:6889-6893 (1979); Kunkel, T.A., In: 
Nucleic Acids and Molecular Bioloav . Vol. 2, 124-135 (Eckstein, F. 
et al .. eds.), Springer-Verlag, Berlin, (1988); Olsen, D.B. et al .. Proc. 
Natl. Acad. Sci. ru. S.A.I 87:1451-1455 (1990); Griep, M.A. et al .. 
Biochemistrv 29 :9006-9014 (1990); Sayers, J.R. et al .. Nucl. Acids 
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Res. 1 6:791-802 (1988)). Alternatively, phosphorothioate 
nucleotide derivatives can be Incorporated synthetically into an 
oligonucleotide (Zon, G. et al .. Anti-Canc. Drug Pes. 6.:539-568 
(1991)). 

The primer molecules are permitted to hybridize to a 
complementary target nucleic acid molecule, and are then extended, 
preferably via a polymerase, to form an extension product. The 
presence of the phosphorothioate nucleotides in the primers 
renders the extension product resistant to nuclease attack. As 
indicated, the amplification products containing phosphorothioate 
or other suitable nucleotide derivatives are substantially resistant 
to "elimination" (i.e. degradation) by "5'-»3'" exdnucleases such as 
T7 exonuclease or exonuclease, and thus a 5'-^ 3' exonuclease will 
be substantially incapable of further degrading a nucleic acid 
molecule once it has encountered a phosphorothioate residue. 

Since the target molecule lacks nuclease resistant residues, 
the incubation of the extension product and its template - the 
target - in the presence of a 5'-»3' exonuclease results in the 
destruction of the template strand, and thereby achieves the 
preferential production of the desired single strand. 

D. Solid Phase Attachment of DNA 

The preferred method of determining the identity of the 
polymorphic site of a polymorphism involves nucleic acid 
hybridization. Although such hybridization can be performed in 
solution (Berk, A.J., et al. £eJi 12 :721 -732 (1977); Hood, LE., et al. . 
In: Molecular Bioloav of Eukarvotic Cells: A Problems Approach . 
Menio Park, CA: Benjamin-Cummings, (1975); Wetmer, J.G., 
Hybridization and Renaturation Kinetics of Nucleic Acids . Ann. Rev. 
Biophvs. Bioena. 5:337-361 (1976); Itakura, K., et al.. Ann. Rev. 
Biochem. 53 :323-356. (1984)), it is preferable to employ a solid- 
phase hybridization assay (see, Saiki, R.K. et al. . Proc. Natl. Acad. 
Sci. (U.S.A.) 86:6230-6234 (1989); Gilham et al. . J. Amer. Chem. 
Soc. 86 :4982 (1964) and Kremsky et al. . Nucl. Acids Res. 
15:3131-3139 (1987)). 
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Any of a variety of methods can be used to immobilize 
oligonucleotides to the solid support. One of the most widely used 
methods to achieve such an immobilization of oligonucleotide 
primers for subsequent use in hybridization-based assays consists 
of the non-covaient coating of these solid phases with streptavidin 
or avidin and the subsequent immobilization of biotinylated 
oligonucleotides (Holmstrom, K. et al.. Anal. Biochem. 209 :278-2fl.'^ 
(1993)). Another known method (Running. J.A. et al.. BioTechnig ues 
8:276-277 (1990); Newton, C.R. et al. Nucl. Acids Res. 21:1M^Pi. 
1162 (1993)) requires the pre-coating of the polystyrene or glass 
solid phases with poly-L-Lys or poly L-Lys, Phe, followed by the 
covalent attachment of either amino- or sulfhydryl-modified 
oligonucleotides using bifunctional crosslinking reagents. Both 
methods have the disadvantage of requiring the use of modified 
oligonucleotides as well as a pre-treatment of the solid phase. 

In another published method (Kawai, S et al. . Anal. Biochem. 
209.-63-69 (1993)), short , oligonucleotide probes were ligated 
together to form multimers and these were ligated into a phagemid 
vector. Following in vitro amplification and isolation of the 
single-stranded form of these phagemids, they were immobilized 
onto polystyrene plates and fixed by UV irradiation at 254 nm. The 
probes immobilized in this way were then used to capture and 
detect a biotinylated PGR product. 

A method for the direct covalent attachment of short, 5'- 
phosphorylated primers to chemically modified polystyrene plates 
("Covalink" plates, Nunc) has also been published (Rasmussen, S.R. 
et al. . Anal. Biochem. 1 98 :138-142 (1991)). The covalent bond 
between the modified oligonucleotide and the solid phase surface 
is introduced by condensation with a water-soluble carbpdiimide. 
This method is claimed to assure a predominantly 5'-attachment of 
the oligonucleotides via their 5'-phosphates; however, it requires 
the use of specially prepared, expensive plates. 

Most preferably, such immobilization of oligonucleotides 
(preferably between 15 and 30 bases) is accomplished using a 
method that can be used directly, without the need for any pre- 
treatment of commercially available polystyrene microwell plates 
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(ELISA plates) or microscope glass slides. Since 96 well 
polystyrene plates are widely used in ELISA tests, there has been 
significant Interest In the development of methods for the 
immobilization of short oligonucleotide primers to the wells of 
these plates for subsequent hybridization assays. Also of interest 
is a method for the immobilization to microscope glass slides, 
since the latter are used in the so-called Slide Immunoenzymatic 
Assay (SIA) (de Macario, E.G. et al. . BioTechniguRR 
(1985)). 

The solid support can be glass, plastic, paper, etc. The 
support can be fashioned as a bead, dipstick, test tube, etc. In a 
preferred embodiment, the support will be a microtiter dish, having 
a multiplicity of wells. The conventional 96-well microtiter 
dishes used in diagnostic laboratories and in tissue culture are a 
preferred support. The use of such a support allows the 
simultaneous determination of a large number of samples and 
controls, and thus facilitates the analysis. Automated delivery 
systems can be used to provide reagents to such microtiter dishes. 
Similarly, spectrophotometric methods can be used to analyze the 
polymorphic sites, and such analysis can be conducted using 
automated spectrophotometers. 

One aspect of the present invention concerns a method for 
Immobilizing oligonucleotides for such analysis. In accordance 
with the method, any of a number of commercially available 
polystyrene plates can be used directly for the immobilization, 
provided that they have a hydrophilic surface. Examples of suitable 
plates include the Immulon 4 plates (Dynatech) and the Maxisorp 
plates (Nunc). The Immobilization of the oligonucleotides to the 
plates is achieved simply by incubation in the presence of a 
suitable salt. No immobilization takes place in the absence of a 
salt, i.e., when the oligonucleotide is present in a water solution. 
Examples for suitable salts are: 50-250 mM NaCI; 30-100 mM 1- 
ethyl-3-(3'-dimethylamlnopropyl)carbodiimide hydrochloride (EDC), 
pH 6.8; 50-150 mM octyldimethylamine hydrochloride, pH 7.0; 50- 
250 mM tetramethylammonium chloride. The immobilization is 
achieved by Incubation, preferably at room temperature for 3 to 24 
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hours. After such incubation, the plates are washed, preferably 
with a solution Q^JO mM.Tris HCI, pH 7.5, contaimrig 150 mM NaCI 
and 0.05% vol.^S»d*^rff^ffiM^^ 

important role of blocking all free oligonucleotide binding sites 
still present on the polystyrene surface, so that no nonspecific 
binding of oligonucleotides can take place during the subsequent 
hybridization steps. Using radioactively labeled oligonucleotides, 
the amount of immobilized oligonucleotides per well was 
determined to be at least 500 fmoles. The oligonucleotides are 
immobilized to the surface of the plate with sufficient stability 
and can only be removed by prolonged incubations with 0.5 M NaOH 
solutions at elevated temperatures. No oligonucleotide is removed 
by washing the plate with water, TNTw (Tween 20), PBS, 1.5 M 
NaCI, or other similar solutions. 

The immobilized oligonucleotides can be used to capture 
specific DNA sequences by hybridization. The hybridization is 
usually carried out in a solution containing 1.5 M .NaCI and 10 mM 
EDTA, for 15 to 30 minutes at room temperature. Other 
hybridization conditions can also be used. More than 400 fmoles of 
a specific DNA sequence was found to hybridize to the immobilized 
oligonucleotide in one well. This DNA is bound to the initially 
immobilized oligonucleotide only via Watson-Crick hydrogen bonds 
can be easily removed from the wells by a brief wash with a 0.1 M 
NaOH solution, without removing the initially attached 
oligonucleotide from the plate. If the captured DNA fragment is 
nonradioactively labeled, e.g., with a biotin residue, the detection 
can be carried out using a suitable enzyme-linked assay. 

Although no modifications have to be introduced into the 
synthetic oligonucleotides, the method also allows for the 
immobilization of labeled (e.g., biotinylated) oligonucleotides, if 
desired. The amount of oligonucleotide that can be immobilized In 
a single well of an ELISA plate by this method is at least 500 
fmoles. The oligonucleotides thus immobilized onto the solid phase 
can hybridize to suitable templates and also participate in 
enzymatic reactions like template-directed extensions and 
ligations. 
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For high volume testing applications, it is desirable to use 
non-radioactive detection methods. Thus, the use of haptenated 
dideoxynucleotides is preferred; the use of biotinylated 
dideoxynucleotides is particularly preferred as such modification 
would render the incorporated base detectable by the standard 
avidin (or streptavidin) enzyme conjugates used in ELISA assays. 
The biotinylated ddNTPs are preferably prepared by reacting the 
four respective (3-aminopropyn-1 -yl)nucleoside triphosphates 
with sulfosuccinimidyl 6-(biotinamido)hexanoate. Thus, (3- 
aminopropyn-1-yl) nucleoside 5'-triphosphates are prepared as 
described by Hobbs, F.W. (iL_Qrs^_CherrL M:3420-3422 (1989)) and 
by Hobbs, F.W. et al . (U.S. Patent No. 5,047,519). The (3- 
aminopropyn-1-yl)nucleoside 5'-triphosphate (50 mol) is dissolved 
in 1 ml of pH 7.6, 1 M aqueous triethylammonium bicarbonate 
(TEAB). Sulfosuccinimidyl 6-(biotinamido) hexanoate sodium salt 
(Pierce, 55.7 mg, 100 mol) is added and the solution is heated to 
50°C in a stoppered tube for 2 hr. The reaction mixture is diluted 
to 10 ml with water and applied to a DEAE-Sephadex A-25-120 
column (1.6 x 19 cm). The column is eluted with a linear gradient 
of pH 7.6 aqueous TEAB (0.1 M to 1.0 M) and the eluent monitored at 
270 nm. The late-eluting major peak is collected, stripped, and co- 
evaporated with ethanol. The crude product, containing 
biotinylated nucleoside triphosphate and, in some cases, 
contaminating starting material, is further purified by reverse 
phase column chromatography (Baker C-18 packing, 2 x 12 cm bed). 
The material is loaded in 0.1 M pH 7.6 TEAB and eluted with a step 
gradient of acetonitriie in 0.1 M pH 7.6 TEAB (0% to 36%, 2% 
increments, 8 ml/step). In all cases, the biotinylated product Is 
more strongly retained and cleanly resolved from the starting 
material. Product-containing fractions are pooled, stripped, and 
co-evaporated with ethanol. The product is taken up in water and 
the yield calculated using the absorption coefficient for the 
starting nucleotide. The 3h NMR and 31 p nMR spectra are 
consistent with the expected structure and confirm the absence of 
phosphorus containing or nucleotide-derived impurities. The 
materials are observed to be >99% pure by HPLC (Waters Bondapak 
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C-18, 4.6 X 250 mm, 1 ml/min, 1 to 35% CHgCN/pH 7/0.01 M 
triethylammonium acetate). 

The synthesis of 5-{3-(6-biotinamido(hexanoamido) propyn- 
1 -yl)-2',3'-dideoxyuridine-5'-triphosphate has an approximate 
yield of 25% (assuming = 12,400 at 291.5 nm); HPLC tx = 16.1 
min. 

The synthesis of 5-{3-(6-biotinamido(hexanoamido) propyn- 
1-yl)-2',3'-dideoxycytidine-5'-triphosphate has an approximate 
yield of 63% (assuming = 9,230 at 294.5 nm); HPLC tx = 19.4 min. 

The synthesis of 7-(3-(6-biotinamido(hexanoamido) propyn- 
1-yl)-7-deaza-2',3'-dideoxyadenosine-5'-triphosphate has an 
approximate yield of 39% (assuming = 13,600 at 278.5 nm); HPLC 
tx = 23.1 min. 

The synthesis of 7-(3-(6-biotinamido(hexanoamido) propyn- 
1-yl)-7-deaza-2',3'-dideoxyguanosine-5'-triphosphate has an 
approximate yield of 44% (assuming = 9,300 at 291 nm); HPLC tx = 
21.2 min. 



E. Solid Phase Analysis of Polymorphic Sites 

1. Polymerase-Mediated Analysis 

Although the identity of the nucleotide(s) of the polymorphic 
sites of the present invention can be determined in a variety of 
ways, an especially preferred method exploits the oligonucleotide- 
based diagnostic assay of nucleic acid sequence variation disclosed 
by Goelet, P. et al. (PCT Application W092/15712, -hetetn- 
i noorporat o d — ty — ref e r e noo). In this assay, a purified 
oligonucleotide having a defined sequence (complementary to an 
immediate proximal or distal sequence of a polymorphism) is bound 
to a solid support, especially a microtiter dish. A sample, 
suspected to contain the target molecule, or an amplification 
product thereof, is placed in contact with the support, and any 
target molecules present are permitted to hybridize to the bound 
oligonucleotide. 
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In one preferred embodiment, an oligonucleotide having a 
sequence that is complementary to an immediately distal sequence 
of a polymorphism is prepared using the above-described methods 
(and preferably that of Nikiforov, T. (U.S. Patent Application Serial 
No. 08/005,061). The terminus of the oligonucleotide is attached 
to the solid support, as described, for example by Goelet, P. et al 
(PCT Application WO 92/15712), such that the 3'-end of the 
oligonucleotide can serve as a substrate for primer extension. 

The immobilized primer is then incubated in the presence of a 
DNA molecule (preferably a genomic DNA molecule) having a single 
nucleotide polymorphism whose immediately 3'-distal sequence is 
complementary to that of the immobilized primer. Preferably, such 
incubation occurs in the complete absence of any dNTP (i.e. dATP, 
dCTP, dGTP, or dTTP), but only in the presence of one or more chairi 
terminating nucleotide triphosphate derivatives (such as a dideoxy 
derivative), and under conditions sufficient to permit the 
incorporation of such a derivative on to the 3'-terminus of the 
primer. As will be appreciated, where the polymorphic site is such 
that only two or three alleles exist (such that only two or three 
species of dNTPs, respectively, could be incorporated into the 
primer extension product), the presence of unusable nucleotide 
triphosphate(s) in the reaction is immaterial. In consequence of 
the incubation, and the use of only chain terminating nucleotide 
derivatives, a single dideoxynucleotide is added to the 3'-terminus 
of the primer. The identity of that added nucleotide is determined 
by, and is complementary to, the nucleotide of the polymorphic site 
of the polymorphism. 

In this embodiment, the nucleotide of the polymorphic site is 
thus determined by assaying which of the set of labeled 
nucleotides has been incorporated onto the 3'-terminus of the 
bound oligonucleotide by a primer-dependent polymerase. Most 
preferably, where multiple dideoxynucleotide derivatives are 
simultaneously employed, different labels will be used to permit 
the differential determination of the identity of the incorporated 
dideoxynucleotide derivative. 
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2. Polymerase/Ligase-Mediated Analysis 

In an alternative embodiment, the identity of the nucleotide 
of the polymorphic site is determined using a polymerase/ligase- 
mediated process. As in the above embodiment, an oligonucleotide 
primer is employed, that is complementary to the immediately 3'- 
distal invariant sequence of the SNP. A second oligonucleotide, is 
tethered to the solid phase via its 3'-end. The sequence of this 
oligonucleotide is complementary to the 5'-proximal sequence of 
the polymorphism being analyzed, but is incapable of hybridizing to 
the oligonucleotide primer. 

These oligonucleotides are incubated in the presence of DNA 
containing the single nucleotide polymorphism that is to be 
analyzed, and at least one 2', 5'-deoxynucleotide triphosphate. The 
incubation reaction further includes a DNA polymerase and a DNA 
ligase. Thus, for example, where the polymorphism of clone 177-2 
(Table 1 ) is being evaluated, and the tethered oligonucleotide' could 
comprise the 3'-distal sequence of SEQ ID N0:2, the second 
oligonucleotide would have the 5'-proximal sequence of SEQ ID 
N0:1. 

The tethered and soluble oligonucleotides are thus capable of 
hybridizing to the same strand of the single nucleotide 
polymorphism under analysis. The sequence considerations cause 
the two oligonucleotides to hybridize to the proximal and distal 
sequences of the SNP that flank the polymorphic site (X) of the 
polymorphism; the hybridized oligonucleotides are thus separated 
by a "gap" of a single nucleotide at the precise position of the 
polymorphic site. 

The presence of a polymerase and a 2', 5'-deoxynucleotide 
triphosphate complementary to (X) permits ligation of the primer 
extended with the complementary 2', 5'-deoxynucleotide 
triphosphate to the immobilized oligo complementary to the distal 
sequence, a 2', 5'-deoxynuc!eotide triphosphate that is 
complementary to the nucleotide of the polymorphic site permits 
the creation of a ligatable substrate. The ligation reaction 
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immobilizes the 2', S'-deoxynucieotide and the previously soluble 
primer oligonucleotide to the solid support. 

The identity of the polymorphic site that was opposite the 
"gap" can then be determined by any of several means. In a 
preferred embodiment, the 2', 5'-deoxynucleotide triphosphate of 
the reaction is labeled, and its detection thus reveals the identity 
of the complementary nucleotide of the polymorphic site. Several 
different 2', 5'-deoxynucleotide triphosphates may be present, each 
differentially labeled. Alternatively, separate reactions can be 
conducted, each with a different 2', 5'-deoxynucieotide 
triphosphate. In an alternative sub-embodiment, the 2', 5'- 
deoxynucleotide triphosphates are unlabeled, and the second, 
soluble oligonucleotide is labeled. Separate reactions are 
conducted, each using a different unlabeled 2', 5'-deoxynucleotide 
triphosphate. The reaction that contains the complementary 
nucleotide permits the ligatable substrate to form, and is detected 
by detecting the immobilization of the previously soluble 
oligonucleotide. 

F. Signal-Amplification 

The sensitivity of nucleic acid hybridization detection assays 
may be increased by altering the manner in which detection is 
reported or signaled to the observer. Thus, for example, assay 
sensitivity can be increased through the use of detectably labeled 
reagents. A wide variety of such signal amplification methods 
have been designed for this purpose. Kourilsky et al. (U.S. Patent 
4,581,333) describe the use of enzyme labels to increase 
sensitivity in a detection assay. Fluorescent labels (Albarella et 
aL, EP 144914), chemical labels (Sheldon III et al. . U.S. Patent 
4,582,789; Albarella et al. . U.S. Patent 4,563,417), modified bases 
(Miyoshi et al., EP 119448), etc. have also been used in an effort to 
improve the efficiency with which hybridization can be observed. 

It is preferable to employ fluorescent, and more preferably 
chromogenic (especially enzyme) labels, such that the identity of 
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the incorporated nucleotide can be determined in an automated, or 
semi-automated manner using a spectrophotometer. 

IV. The Use of SNP Genotyping in Methods of Genetic 
Analysis 

A. General Considerations for Using Single 
Nucleotide Polymorphisms in Genetic Analysis 

The utility of the polymorphic sites of the present invention 
stems from the ability to use such sites to predict the statistical 
probability that two individuals will have the same alleles for any 
given polymorphisms. 

Statistical analysis of SNPs can be used for any of a variety 
of purposes. Where a particular animal has been previously tested, 
such testing can be used as a "fingerprint" with which to determine 
if a certain animal is, or is not that particular animal. 

Where a putative parent or both parents of an individual have 
been tested, the methods of the present invention may be used to 
determine the likelihood that a particular animal is or is not the 
progeny of such parent or parents. Thus, the detection and analysis 
of SNVs can be used to exclude patemity of a male for a particular 
individual (such as a stallion's paternity of a particular foal), or to 
assess the probability that a particular individual is the progeny of 
a selected female (such as a particular foal and a selected mare). 

As indicated below, the present invention permits the 
construction of a genetic map of a target species. Thus, the 
particular array of polymorphisms identified by the methods of the 
present invention can be correlated with a particular trait, in order 
to predict the predisposition of a particular animal (or plant) to 
such genetic disease, condition, or trait. As used herein, the term 
"trait" is intended to encompass "genetic disease," "condition," or 
"characteristics." The term, "genetic disease" denotes a 
pathological state caused by a mutation, regardless of whether 
that state can be detected or is asymptomatic. A "condition" 
denotes a predisposition to a characteristic (such as asthma, weak 




bones, blindness, ulcers, cancers, heart or cardiovascular illnesses, 
skeleto-muscular defects, etc.). A "characteristic" is an attribute 
that imparts economic value to a plant or animal. Examples of 
characteristics include longevity, speed, endurance, rate of aging, 
5 fertility, etc. 

B. Identification and Parentage Verification 

The most useful measurements for determining the power of 
1 0 an identification and paternity testing system are: (i) the 
"probability of identity" (p(ID)) and (ii) the "probability of 
exclusion" (p(exc)). The p(ID) calculates the likelihood that two 
random individuals will have the same genotype with respect to a 
given polymorphic marker. The p(exc) calculates the likelihood, 

1 5 with respect to a given polymorphic marker, that a random male 

will have a genotype incompatible with him being the father in an 
average paternity case in which the identity of the mother is not in 
question. Since single genetic loci, including loci with numerous 
alleles such as the major histocompatibility region, rarely provide 

2 0 tests with adequate statistical confidence for paternity testing, a 

desirable test will preferably measure multiple unlinked loci in 
parallel. Cumulative probabilities of identity or non-identity, and 
cumulative probabilities of paternity exclusion are determined for 
these multi-locus tests by multiplying the probabilities provided 

2 5 by each locus. 

The statistical measurements of greatest interest are: (i) the 
cumulative probability of non-identity (cum p(nonlD)), and (ii) the 
cumulative probability of paternity exclusion (cum p(exc)). 

The formulas used for calculating these probability values 
30 are given below. For simplicity these are given first for 2-allele 
loci, where one allele is termed type A and the other type B. In 
such a model, four genotypes are possible: AA, AB, BA, and BB 
(types AB and BA being indistinguishable biochemically). The 
allelic frequency is given by the number of times A (f(A), the 

3 5 frequency of A is denoted by "p") or B (f(B), the frequency of B is 
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denoted by "q," where q = 1-p) is found in tlie haploid genome. The 
probability of a given genotype at a given locus: 

Homozygote: p{AA)= 

Single Heterozygote: p{AB)= p{BA)= pq = p(1-p) 
Both Heterozygotes: p{AB+BA)= 2pq = 2p(1-p) 

Homozygote: p[BB)= <f - (l-p)^ 

The probability of identity at one locus (i.e the probability 
that two individuals, picked at random from a population will have 
identical genotypes at a given locus) is given by the equation: 

p(/D) = (p2)2 + (2pc/)2 + (q2)2 

The cumulative probability of identity for n loci is therefore 
given by the equation: 

cum p{ID) = ^p{ID,)p{ID^)p{ID,)....p{ID„) 

The cumulative probability of non-identity for n loci (i.e. the 
probability that two individuals will be different at 1 or more loci) 
is given by the equation: 

cum p{nonlD) = 1 - cum p{!D) 
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The probability of parentage exclusion (representing the 
probability that a random nnale will have a genotype, with respect 
to a given iocus, that nnakes him incompatible as the sire in an 
average paternity case where the identity of the mother is not in 
question) is given by the equation: 

p{exc) = pqO'pq) 

The probability of non-exclusion (representing the probability 
at a given locus that a random male will not be biochemically 
excluded as the sire in an average paternity case) is given by the 
equation: 

pinon-exc) = 1 - p{exc) 

The cumulative probability of non-exclusion (representing the 
value obtained when n loci are used) is thus: 

cum p{non'exc) = ^p{non'exc^)p{non'exc^)p{non-exc^)....p{non'Bxc^) 



The cumulative probability of exclusion (representing the 
probability, using a panel of n loci, that a random male will be 
biochemically excluded as the sire in an average paternity case 
where the mother is not in question) is given by the equation: 

cum p{exc) = 1 - cum pinon-exc) 

These calculations may be extended for any number of alleles 
at a given locus. For example, the probability of identity p(lD) for a 
3-allele system where the alleles have the frequencies in the 




population of p, q and r, respectively, is equal to the sum of the 
squares of the genotype frequencies: 

p{ID) = p4 + (2pq)2 + {2qr)^ + (2pr)2 + + 

Similarly^ the probability of exclusion for a three allele 
system is given by: 

p{exc) = pqO'Pq) + qr^-qr) + pr^-pr) + Spqr^-pqr) 

In a locus of n alleles, the appropriate binomial expansion is 
used to calculate p(ID) and p(exc). 

Figures 4 and 5 show how the cum p(nonlD) and the 
cum p(exc) increase with both the number and type of genetic loci 
used. It can be seen that greater discriminatory power is achieved 
with fewer markers when using three allele systems. In Figures 4 
and 5, the triangles trace the increase in probability values with 
increasing numbers of loci with two alleles where the common 
allele is present at a frequency of p = 0.79, The crosses in Figures 
4 and 5 show the same analysis for increasing numbers of three- 
allele loci where p = 0,51, q = 0.34 and r = 0.15. 

The choice between whether to use loci with 2, 3 or more 
alleles is however largely influenced by the above-described 
biochemical considerations. A polymorphic analysis test may be 
designed to score for any number of alleles at a given locus. If 
allelic scoring is to be performed using gel electrophoresis, each 
allele should be easily resolvable by gel electrophoresis. Since the 
length variations in multiple allelic families are often small, 
human DNA tests using multiple allelic families include statistical 
corrections for mistaken identification of alleles. Furthermore, 
although the appearance of a rare allele from a multiple allelic 
system may be highly informative, the rarity of these alleles 
makes accurate measurements of their frequency in the population 
extremely difficult. To correct for errors in these frequency 
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estimates when using rare alleles, the statistical analysis of this 
data must include a measure of the cumulative effects of 
uncertainty in these frequency estimates. The use of these 
multiple allelic systems also increases the likelihood that new or 
rare alleles in the population will be discovered during the course 
of large population screening. The integrity of previously collected 
genetic data would be empirically revised to reflect the discovery 
of a new allele. 

In view of these considerations, although the use of loci with 
many alleles could potentially offer some short-term advantages 
(because fewer loci would need to be screened), it is preferable to 
perform polymorphic analyses using loci with fewer alleles that 
are: (i) more frequently represented, and (ii) easier to measure 
unambiguously. Tests of this type can achieve the same power of 
discrimination as tests based on more highly polymorphic loci, 
provided the same total number of alleles is collected from a 
series of unlinked loci. 

C. Gene Mapping and Genetic Trait Analysis Using 
SNPs 

The polymorphisms detected in a set of individuals of the 
same species (such as humans, horses, etc.), or of closely related 
species, can be analyzed to determine whether the presence or 
absence of a particular polymorphism correlates with a particular 
trait. 

To perform such polymorphic analysis, the presence or 
absence of a set of polymorphisms (i.e. a "polymorphic array") is 
determined for a set of the individuals, some of which exhibit a 
particular trait, and some of which exhibit a mutually exclusive 
characteristic (for example, with respect to horses, brittle bones 
vs. non-brittle bones; maturity onset blindness vs. no blindness; 
predisposition to asthma, cardiovascular disease vs. no such 
predisposition). The alleles of each polymorphism of the set are 
then reviewed to determine whether the presence or absence of a 
particular allele is associated with the particular trait of interest. 
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Any such correlation defines a genetic map of tine individual's 
species. Alleles that do not segregate randomly with respect to a 
trait can be used to predict the probability that a particular animal 
will express that characteristic. For example, if a particular 
5 polymorphic allele is present in only 20% of the members of a 
species that exhibit a cardiovascular condition, then a particular 
member of that species containing that allele would have a 20% 
probability of exhibiting such a cardiovascular condition. As 
indicated, the predictive power of the analysis Is increased by the 

1 0 extent of linkage between a particular polymorphic allele and a 
particular characteristic. Similarly, the predictive power of the 
analysis can be increased by simultaneously analyzing the alleles 
of multiple polymorphic loci and a particular trait. In the above 
example, if a second polymorphic allele was found to also be 

o 1 5 present in 20% of members exhibiting the cardiovascular condition, 
however, all of the evaluated members that exhibited such a 
cardiovascular condition had a particular combination of alleles for 
these first and second polymorphisms, then a particular member 
containing both such alleles would have a very high probability of 

2 0 exhibiting the cardiovascular condition. 
The detection of multiple polymorphic sites permits one to 

define the frequency with which such sites independently 
segregate in a population. If, for example, two polymorphic sites 
p segregate randomly, then they are either on separate chromosomes, 

2 5 or are distant to one another on the same chromosome. Conversely, 
two polymorphic sites that are co-inherited at significant 
frequency are linked to one another on the same chromosome. An 
analysis of the frequency of segregation thus permits the 
establishment of a genetic map of markers. Thus, the present 
30 invention provides a means for mapping the genomes of plants and 
animals. 

The resolution of a genetic map is proportional to the number 
of markers that it contains. Since the methods of the present 
invention can be used to isolate a large number of polymorphic 
35 sites, they can be used to create a map having any desired degree of 
resolution. 
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The sequencing of the polymorphic sites greatly increases 
their utility in gene nnapping. Such sequences can be used to design 
oligonucleotide primers and probes that can be employed to "walk" 
down the chromosome and thereby identify new marker sites 
(Bender, W. et al. , J. Supra. Molec, Struc. 1 0fsuppl,k 32 (1979); 
Chinault, A.C. et ah . Gene 5:111-126 (1979); Clarke, L. et al, . Nature 
287:504-509 (1980)), 

The resolution of the map can be further increased by 
combining polymorphic analyses with data on the phenotype of 
other attributes of the plant or animal whose genome is being 
mapped. Thus, if a particular polymorphism segregates with brown 
hair color, then that polymorphism maps to a locus near the gene or 
genes that are responsible for hair color. Similarly, biochemical 
data can be used to increase the resolution of the genetic map. In 
this embodiment, a biochemical determination (such as a serotype, 
isoform, etc.) is studied in order to determine whether it co- 
segregates with any polymorphic site. Such maps can be used to 
identify new gene sequences, to identify the causal mutations of 
disease, for example. 

Indeed, the identification of the SNPs of the present 
invention permits one to use complimentary oligonucleotides as 
primers in PGR or other reactions to isolate and sequence novel 
gene sequences located on either side of the SNP. The invention 
includes such novel gene sequences. The genomic sequences that 
can be cionaliy isolated through the use of such primers can be 
transcribed into RNA, and expressed as protein. The present 
invention also includes such protein, as well as antibodies and 
other binding molecules capable of binding to such protein. 

The invention is illustrated below with respect to two of its 
embodiments - horses and humans. However, because the 
fundamental tenets of genetics apply irrespective of species, such 
illustration is equally applicable to any other species. Those of 
ordinary skill would therefore need only to directly employ the 
methods of the above invention to isolate SNPs in any other 
species, and to thereby conduct the genetic analysis of the present 
invention. 




As indicated above, LOD scoring methodology has been 
developed to permit the use of RFLPs to both track the inheritance 
of genetic traits, and to construct a genetic map of a species 
(Lander, S. et al. , Proc. Natl. Acad. Sci. fU.S.A,^ 83:7353-7.^^7 

(1986) ; Lander, S. et a!. . Proc. Natl. Acad. Sci, (U.S.A.) 84:2363-2367 

(1987) ; Donis-Keiler, H. et al. . Cell 51:319-337 (1987); Lander, 8, 
et aL . Genetics 121 :185-199 (1989)). Such methods can be readily 
adapted to permit their use with the polymorphisms of the present 
invention. Indeed, such polymorphisms are superior to RFLPs and 
STRs in this regard. Due to the frequency of SNPs, it is possible to 
readily generate a dense genetic map. Moreover, as indicated 
above, the polymorphisms of the present invention are more stable 
than typical (VNTR-type) RFLP polymorphisms. 

The polymorphisms of the present invention comprise direct 
genomic sequence information and can therefore be typed by a 
number of methods. In an RFLP or STR-dependent map, the analysis 
must be gel-based, and entail obtaining an electrophoretic profile 
of the DNA of the target animal. In contrast, an analysis of the 
polymorphisms (SNPs) of the present invention may be performed 
using spectrophotometric methods, and can readily be automated to 
facilitate the analysis of large numbers of target animals. 

Having now generally described the invention, the same will 
be more readily understood through reference to the following 
examples of the isolation and analysis of equine polymorphisms 
which are provided by way of illustration, and are not intended to 
be limiting of the present invention. 

EXAMPLE 1 

DISCOVERY OF EQUINE POLYMORPHISMS 

As an initial step in the identification of equine 
polymorphisms, small shotgun libraries were prepared from 
genomic DNA isolated from peripheral blood leukocytes which had 
been purified on a Ficoll-hypaque density gradient from the blood 
of a single, 15 year old thoroughbred gelding (John Henry). This 
DNA was simultaneously digested to completion with Bam HI and 




Pst i and either used directly or after size fractionation on agarose 
gels. 

Vector pLT14 (a variant of the Stratagene plasmid pKSM13(- 
)) was digested with Bam HI and Pst I and linearized DNA was 
purified from an' agarose gel. For both vector and size-fractionated 
genomic DNA, agarose plugs were solubilized in saturated sodium 
iodide and the DNA was subsequently immobilized on glass powder. 
After washing, the DNA was eluted with water and ethanol 
precipitated with glycogen carrier. 

Ligations with varying vector/insert ratios were effectuated 
with T4 DNA ligase at 4°C. E. coli strain XLI was transformed with 
ligation mixtures and plated on LB agar containing 100 g/ml 
ampicillin. Approximately 50,000 clones were generated in several 
different experiments using size fractionated or unfractionated 
insert DNA. Unplated transformed ceils were stored at -70°C in 7% 
DMSO. Colonies were streaked for isolation and small scale 
plasmid preparations were performed to determine the size of 
inserted equine DNA. Larger scale preparations were performed 
with Qiagen chromatography. 

The sequence of the first 200-300 nucleotides of the genomic 
insert was determined by the chain terminating dideoxynucleoside 
method with T7 DNA polymerase from primers complementary to 
plasmid sequences. This information was used to design synthetic 
oligonucleotide primers complementary to the equine sequence to 
be employed in PGR reactions. 

In most cases, two sets of PGR primers (generally 25-mers) 
were synthesized. The first set was used to amplify, under a 
standardized set of conditions, from genomic DNA. The products of 
these reactions were diluted and used as template DNA in a second 
PGR using nested primers slightly internal to the original set. The 
products of these two reactions were compared to those obtained 
using the original plasmid DNA as template, in most cases, It was 
possible to obtain high quality, single-species products using this 
procedure with no attempt to optimize reaction conditions for any 
particular pair of primers. 



Two different methods were used to screen amplified DNA 
from horses for polymorphic sequences. Initially, PGR fragments 
from a panel of 6 horses were digested with a panel of restriction 
endonucleases having 4 base recognition sites. The products of 
5 these reactions were analyzed by acrylamide gel electrophoresis on 
5% - 7.5% non-denaturing gels. Digestion products which showed 
variability when hybridized to different members of the panel were 
subjected to DNA sequence analysis. Later, DNA sequencing was 
used directly to screen for polymorphic sites. The PGR fragments 
1 0 from five unrelated horses were electroeluted from acrylamide 
gels and sequenced using repetitive cycles of thermostable Taq 
polymerase reaction in the presence of a mixture of dNTPs and 
fluorescent ddNTPs, The products were then separated and 
analyzed using the automated DNA sequencing instrument of 
p 1 5 Applied Biosystems, Inc. The data was analyzed using ABl 
yri software. Differences between sequences of different animals 

ji; were identified by the software and confirmed by inspecting the 

gri relevant portion of the chromatograms on the computer screen. 

^ Differences were concluded to be a DNA polymorphism only if the 

y 20 data was available for both strands, and/or present in more than 
one haploid example among the five horses tested. 

y i 

g EXAMPLE 2 

5 CHARACTERIZATION OF EQUINE POLYMORPHISMS 



25 

The program of identification and characterization of 
polymorphic DNA sequences in randomly selected fragments was 
continued such that approximately 550 plasmids have been 
characterized to this level. The sequences adjacent to the cloning 
30 sites was determined for 200 of these plasmids. Inserts of these 
sequenced plasmids ranged in size from 0.25 to 3.5 kb. Using this 
sequence information, oligonucleotide primers were designed to 
enable PCR amplification of the same genomic region from 
different horses. 

35 In order to identify the nucleotides present at polymorphic 

sites, PCR fragments from 5 horses were purified from acrylamide 



- 48 - 



gels by electroelution and completely sequenced using Taq 
polymerase "Cycle" sequencing biochemistry and automated 
sequencing equipment. Results from the 5 horses were analyzed by 
computer and visually confirmed, DNA sequence variants 
discovered by this method were scored only if the sequence was 
obtained on both strands and the variant sequence had been found in 
more than one haploid example. The 18 clones of Table 1 comprise 
a subset of identified SNPs. In Table 1, the immediately 5'- 
proximal sequence, the identity of the nucleotide of the 
polymorphic site, and the immediately 3'-distal sequence of each 
SNP is presented. For each SNP, Such sequences are shown in the 
horizontal rows. The sequences of double-stranded DNA in Table 1 
is presented in compliance with the Sequence Listing requirements 
of the United States Patent and Trademark Office, Thus, all 
sequences are presented in the same orientation (5'->3'). The 
organization of the Table is illustrated in Figure 6 with respect to 
an illustrative SNP, clone 177-2, This SNP has a polymorphic site 
capable of having either a C or a T in one strand, and a G or A in the 
opposite strand. The 5'-proximal DNA sequence that immediately 
precedes the polymorphic site in the C/T strand is designated as 
SEQ ID N0:1. The 3'-distal sequence that immediately follows the 
polymorphic site in the C/T strand is designated as SEQ ID N0:2. 
The 5'-proximal DNA sequence that immediately precedes the 
polymorphic site in the G/A strand is designated as SEQ ID N0:3. 
The 3'-distal sequence that immediately follows the polymorphic 
site in the G/A strand is designated as SEQ ID N0:4. Bearing in 
mind that the sequences are written in the same orientation 
(5'-^3'), it will be seen that the sequences of SEQ ID N0:1 and SEQ 
ID N0:4 are complimentary; similarly, the sequences of SEQ ID NO:2 
and SEQ ID N0:3 are complimentary. The sequences that flank a 
particular polymorphic site are thus obtained by combining the 
proximal sequence of one row with the distal sequence also shown 
in the same row. 
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TABLE 1 II 


POLYMORPHIC LOCI || 


SEQ ID 
NO. 




U3 CO 


O CN 
t-l rH 


rH rH 


00 o 
rH (N 


CN ^ 
CN CN 


U3 <X3 
CN CN 


O CN 

m n 


n m 


00 o 


CN 


U3 00 


O CN 

in in 


^ vo 
in in 


00 o 

in VD 


CN ^ 


U3 00 


O CN 


3 ' DISTAL SEQUENCE 


TGCAGAAATTCTAAGGTGTT 
CCCACAGCACTTAGAGCTGC 


TGAGGGAAAAATGATGATGC 
TAGTGGATCATCCCAGAGCT 


AAAGTCAGAACAATGATTGC 
ATGGCCATCAAATTAGTTTT 


TGAGTTCAGCGAGTGTCAGA 
TCCTGTTCTTGGGCCTTGTG 


GTTCCCTTTTGTTTCAGACC 
AAATAATGTGGTCTTTCCCA 


CCATGGAGTCATAGATAAGT 
CCGGATGCTTCTTACTCATG 


AGCGAGTGTCAGAGTTGTGT 
AACTCAATCCTGTTCTTGGG 


GTCCTACAATTGCCAGGAAG 
AAGGCCCCCCATTTCTTGCT 


TGTGTGTGTGTGTATTTGCT 
ATATATATATATTGATATTC 


GTTTGTTTTAAATTTTGTGA 
AATACAGGCTTAATTATGGC 


ATGTCTTATACCTCCCTTTT 
GTCCTTGAACTACTCTACAC 


CCTGCTGAAGCCTCAGACCG 
AAGGCCTGCTCTCCGTTCAC 


TCAACCTTGCATCATGAGCT 
GGTCATAGTCTAAAGAGCAG 


TCTCCTGCCTTTAGACTCGA 
GACTGAGGTCCCAGCTCAAA 


TTGTTCAGAAGCACAGGTGA 
TATCCACGGCCCAGAGGTTC 


ATCCACTAATGAGGGAAAAA 
ATCCCAGAGCTAGCAAATAC 


AGAGATGCAGCTCTAAGTGC 
CACATCTGTCCCACAACTTC 


GTCCCAATAGTCTGGGATTC 
TTGTGGAGGCTTCCTCATGG 


IDENTIFIED SNP 
ALLELE 




o u 


< £h 


o o 


< Eh 


o o 


Eh < 


o o 


< Eh 


CD CJ 


U C3 


CD O 






CD U 


CD O 


< £h 


O CJ 




< Eh 


U U 




Eh < 


O CJ 


u o 


< Eh 


CD O 


< 


< Eh 


U O 


CD O 




< Eh 


£h < 


CJ CD 


a CD 


5' PROXIMAL SEQUENCE 


GCAGCTCTAAGTGCTGTGGG 
AACACCTTAGAATTTCTGCA 


AGCTCTGGGATGATCCACTA 
GCATCATCATTTTTCCCTCA 


AAAACTAATTTGATGGCCAT 
GCAATCATTGTTCTGACTTT 


CACAAGGCCCAAGAACAGGA 
TCTGACACTCGCTGAACTCA 


TGGGAAAGACCACATTATTT 
GGTCTGAAACAAAAGGGAAC 


CATGAGTAAGAAGCATCCGG 
ACTTATCTATGACTCCATGG 


CCCAAGAACAGGATTGAGTT 
ACACAACTCTGACACTCGCT 


AGCAAGAAATGGGGGGCCTT 
CTTCCTGGCAATTGTAGGAC 


GAATATCAATATATATATAT 
AGCAAATACACACACACACA 


GCCATAATTAAGCCTGTATT 
TCACAAAATTTAAAACAAAC 


GTGTAGAGTAGTTCAAGGAC 
AAAAGGGAGGTATAAGACAT 


GTGAACGGAGAGCAGGCCTT 
CGGTCTGAGGCTTCAGCAGG 


CTGCTCTTTAGACTATGACC 
AGCTCATGATGCAAGGTTGA 


TTTGAGCTGGGACCTCAGTC 1 
TCGAGTCTAAAGGCAGGAGA | 


GAACCTCTGGGCCGTGGATA 11 
TCACCTGTGCTTCTGAACAA || 


GTATTTGCTAGCTCTGGGAT 11 
TTTTTCCCTCATTAGTGGAT || 


GAAGTTGTGGGACAGATGTG 11 
GCACTTAGAGCTGCATCTCT || 


CCATGAGGAAGCCTCCACAA |j 
GAATCCCAGACTATTGGGAC || 




U} ^ Z 








ro in 

tH rH 


rH rH 


rH n 
CN CM 


CN CN 


cry rH 
CN ro 


m in 
m m 


cn 
m m 


rH m 


in 


<7\ rH 

in 


m in 
in in 


in in 


rH fO 


in 


a\ rH 




CLONE 


177-2 


3.-3 


090-2 


324-1 


... 


007-1 




n,-3 




007-3 


... 


085-1 


007-2 


474-1 


178-1 


595-2 


177-1 


CN 

1 

o^ 
in 
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The present specification refers to the above sequences by 
their sequence ID numbers (i.e. SEQ ID NO). To facilitate such 
disclosure, algebraic notation (such as "2n+1") is employed, in 
accordance with conventional algebra. Thus, the designation "SEQ 
5 ID N0:(2n+1)" denotes SEQ ID N0:5 where n=2, and SEQ ID N0:7 where 
n=3, etc. 

EXAMPLE 3 

ALLELIC FREQUENCY ANALYSIS OF EQUINE POLYMORPHISMS IN SMALL 
1 0 POPULATION STUDIES 



Small population studies (50 - 60 animals) of these DNA 
sequence polymorphisms has been carried out on a number of these 
polymorphic sites using Genetic Bit Analysis (GBA), the preferred 
1 5 solid-phase, single nucleotide interrogation system (Goelet, P. et 
aL (WO 92/15712). The 7 steps of the most preferred embodiment 
*le=ntiustrated in Figure 7: 

Step 1: DNA preparation. 

Step 2: Amplification of Target Sequence. After DNA is 
20 prepared from the sample, a specific region of the sample genome 
(locus) is amplified using the PCR. One of the PCR primers is 
modified with four phosphorothioate linkages at the 5'-end. 

Step 3: Exonuclease Digestion and the Generation of Single- 
Stranded Template. The PCR product is digested with exonuclease, 
2 5 leaving the phosphorothioated strand intact. 

Step 4: Hybridization to Capture the Amplified Template. The 
template strand is next hybridized to the appropriate GBA primer 
that is immobilized on the surface of a microtiter well. 

Step 5: Single Base Extension with Polymerase. DNA 
30 polymerase and haptenated ddNTPs are used to extend the GBA 
primer by one base in a template-dependent manner. 

Step 6: Colorimetric detection of the Extension Product. 
After the template is washed away using NaOH, the haptenated base 
is detected using an anti-hapten conjugate and the appropriate 
35 colorimetric substrate. 
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Step 6: Computer-Assisted Interpretation of Genotype. The 
colorimetric data from a number of loci is converted to an SNP 
genotype for the particular individual tested. 

The method is preferably conducted in the following manner: 

GBA Template Preparation. 

Amplification of genomic sequences was performed using the 
polymerase chain reaction (PGR). In a first step, one hundred 
nanograms of genomic DNA was used in a reaction mixture 
containing each first round primer at a concentration of 2 M and 10 
mM Tris pH 8.3, 50 mM KCI, 1.5 mM MgCIg, 0.01% gelatin; and 0.05 
units per I Taq DNA Polymerase (AmpliTaq, Perkin Elmer). 

To obtain single-stranded template for use with solid-phase 
immobilized primer, either of two methods may be used. First, the 
amplification ^ may be mediated using primers that contain 4 



(U.S. patent application serial no. 08/005,061). Alternatively, a 
second round of PGR may be performed using "asymmetric" primer 
concentrations. The products of the first reaction are diluted 
1/1000 in a second reaction. One of the second round primers is 
used at the standard concentration of 2 M while the other is used at 
0.08 M. Under these conditions, single stranded molecules are 
synthesized during the reaction. 

Solid phase immobilization of nucleic acids. 

For the GBA procedure, solid-phase attachment of the 
template-primer complex simplifies washes, buffer exchanges, 
etc., and in principle this attachment can be either via the template 
or the primer. In practice, however, especially when non gel-based 
detection methods are employed, attachment via the primer is 
preferable. This format allows the use of stringent washes (e.g., 
0.2 N NaOH) to remove impurities and reaction side products while 
retaining the haptenated dideoxynucleotide covalently linked to the 
3'-end of the primer. 

Therefore, for GBA reactions in 96-well plates (Nunc Nunclon 
plates, Roskilde, Denmark), the GBA primer was covalently coupled 




derivatives, as taught by Nikiforov, T. 
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to the plate. This was accomplished by incubating 10 pmoles of 
primer having a 5'-amino group per well in 50 of 3 mM sodium 
phosphate buffer, pH 6, 20 mM 1-ethyl-3-(3-dimethylaminopropyl)- 
carbodiimide (EDC) overnight at room temperature. After coupling, 
the plate was washed three times with TNTw. 

GBA in Microwell Plates. 

Hybridization of single-stranded DNA to primers covalently 
coupled to 96-well plates was accomplished by adding an equal 
volume of 3 M NaCI, 20 mM EDTAJo the single-stranded PGR 
product and incubating each well witt^° gO I o f this mixture at 20°C 
for 30 minutes. The plate was subsequently washed three times 
with TNTw. Twenty^vtj of polymerase extension mix containing 
ddNTPs (3 M each, one of which was biotinylated, 5 mM DTT, 7.5 mM 
sodium isocitr^e^^^^__5m^M MnClg, ^94 units perMLc^Klenow DNA 
polymerase) ^ nd ma^c ^d^^/''^mm^^^'yi^[^\S?f^rature. 

Following the extension reaction, the plate was washed once 
with TNTw. Template strands were removed by incubating wells 
with 50 |j.l of 0.2 N NaOH for 5 minutes at room temperature, then 
washing the well with another 50 p.1 of 0.2 N NaOH. The plate was 
then washed three times with TNTw. Incorporation of biotinylated 
ddNTPs was measured by an enzyme-linked assay. Each well was 
incubated with 20 jil of streptavidin-conjugated horseradish 
peroxidase (1/1000 dilution in TNTw of product purchased from 
BRL, Gaithersburg, MD) with agitation for 30 minutes at room 
temperature. After washing 5 times with TNTw, 100 |ii of o- 
phenylenediamine (OPD, 1 mg/ml in 0.1 M citric acid, pH 4.5) (BRL) 
containing 0.012% HgOg was added to each well. The amount of 
bound enzyme was determined kinetically with a Molecular Devices 
model "Vmax" 96-well spectrophotometer. Figures 8A and 88 
illustrate how horse parentage data appears at the microtiter plate 
level. In standard horse parentage testing, samples are arrayed 85 
to a plate (columns 1-11) plus controls (column 12). For each 
horse locus the presence of the two known alleles is determined by 
base specific interrogation on separate plates. The two plates 
shown in figures 8A and 88 are identical in PGR template and GBA 
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primer and differ only in the biotinyiated ddNTP tliat was used In 
the extension reaction (biotin-ddCTP in Figure 8A and biotin-ddTTP 
in Figure 8B). Upon addition of the colorimetric reagent (OPD), the 
absorbance of the resultant color was measured in a Molecular 
Devices microtiter plate reader and the raw data generated in 
milliOD/min per well. The two raw data gray scale representations 
of the absorbance data for these plates are shown in the figures 
arranged in the exact same order as on the microtiter plates. Gray 
scale intensity correlates directly with color production. At this 
biallelic locus the bases detected are C (Figure 8A) and T (Figure 
SB). Approximately 40% of horses tested to date are heterozygotes 
(the sample in well A1, for example) and the remaining homozygous 
for C (A2, for example) or T (83, for example). Synthetic template 
controls include a control C homozygote (well E12), a control T 
homozygots (well F12) and a control heterozygote (well G12). 
Scale refers to milliOD/min at 450 nm. Most positive samples had 
signals above 100 in this case. In this format, for a 28 biallelic 
marker panel horse parentage test, 56 such plates would be 
required for complete typing of the 85 horses. 

Fifty-one random, unrelated horses and three sire/dam/foal 
families were chosen for study in order to establish that a 
reasonable subset of the group of DNA markers found to date was 
likely to provide the desired p(exc) > 0.90, and to assess the power 
of the DNA markers thereby allowing them to be prioritized for 
definitive allelic frequency measurements. 

PGR generated single-stranded template DNA was prepared 
from the genomic DNA of each animal. This material was typed 
with respect to nucleotide variants using GBA. The genotype data 
obtained for each polymorphic site is summarized in Table 2. From 
this genotype data, allelic frequencies were determined and used to 
calculate the p(exc) of each site. The cumulative p(exc) is given 
for the group of 18 sites listed in Tables 1 and 2 is 0.955 for the 
group. In Tables 2-5, the genotype is indicated as either 
homozygote (i.e. PP or QQ) or the heterozygote (PQ). The numbers in 
parentheses denote the number of alleles of the genotype observed. 
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Genotype 3 
Oa (#) 

TT(19) 
TT( 9) 
CC(31) 
GG ( 0) 


GG (17) 
CC (46) 
TT(18) 
GG(ll) 
GG( 9) 
TT( 3) 
GG(5) 
GG(4) 
TT(20) 
GG(9) 
GG (31) 
GG(1) 




Genotype 2 

CT (30) 
CT (24) 
AC (22) 
CG ( 6) 
AT (21) 


I « (N rn M iJi f^^^ ? rn (Ti in CN 
' w Cl. w Ci cilCi £^ (N ro 

O U E-- O O f- O O H O O t5 
< < U < < O < CJ 5 U < < 




Genotype 1 
PP (#) 

CC(ll) 
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AA(35) 
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EXAMPLE 4 

PARENTAGE TESTING 

A family consisting of a sire, dam and offspring was typed 
with respect to the 18 variable sites discussed above with no 
exclusions found. This family had not been previously blood typed. 
Using the preliminary allelic frequency numbers given in Table 2, it 
is possible to construct a p(exc) table pertaining to this specific 
case (Table 3). In general, this Table is constructed assuming that 
the identity of the dam is not in question (although in practice, it 
is possible to exclude the mare if neither of her alleles is inherited 
by the foal). Table 3 shows the typing data for the foal and its dam 
with the sites tested listed in order of informativeness in this 
case. The overall cum p(exc) using 18 loci was 0.942. 











TAB 


LE3 


LOCUS 


FOAL 


DAM 


EXCL'DED 
SIRES 


p(exc) 


p(non- 
exc) 


cum p(non- 
exc) 


cum p(exc) 


459-1 


AC 


CC 


AA 


0.524 


0.476 


0.476 


0.524 


129-1 


AA 


AT 


TT 


0.370 


0.630 


0.300 


0.700 


324-1 


CC 


CT 


TT 


0.321 


0.679 


0.204 


0.796 


595-3 


GG 


GG 


AA 


0.279 


0.721 


0.147 


0.853 


090-2 


GG 


AG 


AA 


0.217 


0.783 


0.115 


0.885 


324-2 


CC 


CT 


TT 


0.151 


0.849 


0.098 


0.902 


595-1 


AA 


AA 


GG 


0.092 


0.818 


0.080 


0.920 


007-3 


AA 


AA 


GG 


0.080 


0.920 


0.073 


0.927 


085-1 


CC 


CC 


GG 


0.071 


0.929 


0.068 


0.932 


474-1 


AA 


AA 


TT 


0.059 


0.941 


0.064 


0.936 


178-1 


AA- 


AG 


GG 


0.043 


0.957 


0.061 


0.939 


595-2 


GG 


GG 


TT 


0.036 


0.964 


0.059 


0.941 


177-1 


CC 


CC 


AA 


0.018 


0.982 


0.058 


0.942 


459-2 


CC 


CC 


GG 


0.003 


0.997 


0.058 


0.942 


007-1 


CG 


CG 




0.000 


1.000 


0.058 


0.942 


007-2 


AG 


AG 




0.000 


1.000 


0.058 


0.942 


177-2 


CT 


CT 




0.000 


1.000 


0.058 


0.942 


177-3 


AG 


AG 




0.000 


1.000 


0.058 


0.942 




EXAMPLE 5 

IDENTITY TESTING 



It is of interest to mal<e use of the population analysis group 
to derive preliminary information concerning other aspects of the 
marker panel. For example, using the allelic frequency data, it is 
possible to calculate a probability of identity [p(ID)] value for the 
18 sites which is equal to 4.79 x 10"' or approximately 1 in 2.1 
million. Thus, one would predict that none of the horses examined 
in the population group would have the same genotype and computer 
analysis of the genotype database revealed this to be the case. As 
shown in Table 4. the p(ID) reaches very small numbers with 
analysis of comparatively few loci. Using the top seven sites, the 
probability of two random animals having different genotypes is 
already 99.9%. 





TABLE 4 


LOCUS 


GENOl'YPH 
1 

PP (#) 


GENOTYPE 
2 

PO (#) 


GENOTYPE 
3 

00 (#) 


P 


q 


P(ID) 


cum 
P(ID) 


177-2 


CC (18) 


CT (23) 


TT (18) 


0.500 


0.500 


0.375 


0.375 


595-3 


AA (14) 


AG (28) 


GG (11) 


0.528 


0.472 


0.376 


0.141 


090-2 


AA (13) 


AG (28) 


GG (17) 


0.466 


0.534 


0.376 


0.053 


324-1 


CC (11) 


CT (30) 


TT (19) 


0.433 


0.567 


0.380 


0.020 


129-1 


AA ( 7) 


AT (33) 


TT (20) 


0.392 


0.608 


0.388 


0.008 


007-1 


AA (22) 


CG (29) 


GG ( 9) 


0.608 


0.392 


0.388 


0.003 


324-2 


CC (21) 


CT (24) 


TT ( 9) 


0.611 


0.389 


0.388 


0.001 


177-3 


AA (26) 


AG (25) 


GG ( 9) 


0.642 


0.358 


0.397 


4.67x10-"^ 


595-1 


AA (25) 


AG (21) 


GG ( 5) 


0.696 


0.304 


0.422 


1.97x10"'^ 


007-3 


AA (27) 


AG (32) 


GG ( 1) 


0.717 


0.283 


0.435 


8.57x10-'^ 


459-1 


AA (5) 


AC (22) 


CC (31) 


0.276 


0.724 


0.440 


3.77x10-5 


085-1 


CC (32) 


CG (24) 


GG ( 4) 


0.733 


, 0.267 


0.447 


1.68x10-5 


007-2 


AA (3) 


AG (25) 


GG (31) 


0.263 


0.737 


0.450 


7.58x10"^ 


474-1 


AA (35) 


AT (21) 


TT ( 4) 


0.758 


0.242 


0.468 


3.55x10-^ 


178-1 


AA (38) 


AG (16) 


GG ( 4) 


0.793 


0.207 


0.505 


1.79x10-6 


595-2 


GG (34) 


GT (13) 


TT ( 3) 


0.810 


0.190 


0,527 


9.45x10-'^ 


177-1 


AA ( 2) 


AC (12) 


CC (46) 


0.133 


0.867 


0.618 


5.84x10-7 


459-2 


CC (53) 


CG ( 6) 


GG ( 0) 


0.949 


0.051 


0-821 4.79x10-7 




False Report Rate 

In the current study, two types of potential false reports can 
be encountered due to either (1) PGR failures or 
(2) incompatibility between the genotype obtained on opposite 
strands. Only data from those animals which had been successfully 
typed in both strands was included in the allelic frequency 
calculations. Sixty horses typed with respect to 18 sites amounts 
to 1,080 genotypings. 95% of all typing experiments were 
successful overall. No typing errors were due to traditional PGR 
failures. 3.8% false reports were encountered at the GBA step 
either because the PGR was unsuccessful at the single strand step 
or due to operator error. 1.1% of all typings produced incompatible 
data between the strands for unknown reasons. 

In sum, the GBA (genetic bit analysis) method is thus a 
simple, convenient, and automatable method for interrogating SNPs. 
In this method, sequence-specific annealing to a solid phase-bound 
primer is used to select a unique polymorphic site in a nucleic acid 
sample, and interrogation of this site is via a highly accurate DNA 
polymerase reaction using a set of novel non-radioactive 
dideoxynucleotide analogs. One of the most attractive features of 
the GBA approach is that, because the actual allelic discrimination 
is carried out by the DNA polymerase, one set of reaction 
conditions can be used to interrogate many different polymorphic 
loci. This feature permits cost reductions in complex DNA tests by 
exploitation of parallel formats and provides for rapid development 
of new tests. 

The intrinsic error rate of the GBA procedure in Its present 
format is believed to be low; the signal-to-noise ratio in terms of 
correct vs. incorrect nucleotide incorporation for homozygotes 
appears to be approximately 20:1. GBA is thus sufficiently 
quantitative to allow the reliable detection of heterozygotes in 
genotyping studies. The presence in the DNA polymerase-mediated 
extension reaction of all four dideoxynucleoside triphosphates as 
the sole nucleotide substrates heightens the fidelity of genotype 
determinations by suppressing misincorporation. GBA can be used 
in any application where point mutation analyses are presently 




employed -- including genetic mapping and linkage studies, genetic 
diagnoses, and identity/paternity testing -- assuming that the 
surrounding DNA sequence is known. 

EXAMPLE 6 

ANALYSIS OF A HUMAN SNP 

Human single nucleotide polymorphisms may be used in the 
same manner as the above-described equine polymorphisms. 
Examples of suitable human polymorphisms are presented in Table 
5. 
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For the purpose of validating the strategy of converting 
human SNPs to a GBA test format, a phenotypically neutral SNP site 
was converted and tested by GBA. This site was selected from the 
Johns Hopkins University 0MB database of human polymorphisms 
The site is met-H on chromosome 7 at q31, mutation position 127 
A to G (Horn, G.T. et al., Clin. Chem. 36, 1614-1619, 1990) The 
following oligonucleotides were synthesized (p=phosphorothioate): 

PGR primer no. 1552 (SEQ ID NO:93) 

5'-CpApTpCpCATGTAGGAGAGCCTTAGTC 

PGR primer no. 1553 (SEQ ID NO:94) 

5'-CCATTTTTGTGTCTTCTAGTCTAAGG 

GBA primer no. 1554 (SEQ ID NO:95) 

5'-TTGAAAGATCGTCAGAAAAATCC 



Human. DNA samples were randomly selected from the DNA 
archives of two families available from the Centre D'Etude du 
Polymorphisme Humaine (CEPH) family collection. A negative 
control, containing no DNA was also used. Sample DNAs were 
amplified by PGR using the above primers and the resulting product 
was analyzed by GBA for two potential bases at the polymorphic 
site, G and A. GBA results were obtained by an endpoint reading of 
absorbance at 450 nm in a microtiter plate reader. The data is 
presented in Table 6. 

Samples 1, 2, 4, 6 and 8 were homozygous for A, samples 7 and 
9 were homozygous for G and samples 3 and 5 were GA 
heterozygotes. These DNAs have not been tested for this biallelism 
by any other method to date. 
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False Rfiport Ratfi 

in the current st^, two types of potential false reports can 
be encountered du\ to either (1) PCR failures or 
(2) mcompatibllity betwee\ the_g^type obtained on opposite 
strands. Only data from Xh^ animai>which had been successfully 

typed m both strands /was\ included \ the allelic frequency 
calcu ations. J^ivtw hr^rSL^ +,,«a^ ...-.^t. \ . . ^ ' 



■ calculations. Sixty hor 
to 1,080 genotypings. 
successful overall. No 
failures. 3.8% false r€ 
either because the RCR 
or due to operator error 



5s typdd with respdct to 18 sites amounts 
95% Vf all typiXg experiments were 
/ping errVs were^jd^re-TTiraditional PCR 
)rts wer^^^Tfcoij^ at the GBA step 
Jnsuccekful at thte single strand step 
VV°/° of a" typings prdduced incompatible 
data between the strands fdAunknown re^ons. \ 

In sum. the GBA (ger^tic bit analyL) meVhod is a simple 
convenient, and automatable n\ethod for inter^gating SNPs In this 
method, sequence-specific annealing to a soli\phase-bound primer 
IS used to select a unique polymorphic site\n a nucleic acid 
sample, and interrogation of this site is via a hi^KHy accurate DNA 
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polyme?ese reaction using a set of novel non-radioactive 
dideoxynOcieotide analogs. One of the most attractive features of 
the GBA ^proach is that, because the actual allelic discrimination 
is carried W by the DNA polymerase, one set of reaction 
conditions ca\i be used to interrogate many different polymorphic 
loci. This feature permits cost reductions in complex DNA tests by 
exploitation of \arallelJflirQ^ and provides for rapid development 
of new tests. ^\ 

The ftntrinsA error rate of Vie GBA procedure in its present 
format is ftelie^ed \o be low; the\signal-to-noise ratio in terms of 
correct vs.\ i/hcorredl nucleotid^.^Worp for homozygotes 

appears toffee appr^xiw*^^ GBA is thus sufficiently 

quantitative Jbjit0T;^rlhe reliable detection of heterozygotes in 
genotyping swSes. TheVresence irl the DNA polymerase-mediated 
extension rea^ion of allVjur dideokynucleoside triphosphates as 
the sole nucleotide substrd^s heightens the fidelity of genotype 
determinations by suppressin^ymisincorporation. GBA can be used 
in any application where poirW mutation analyses are presently 
employed -- including genetic moping and linkage studies, genetic 
diagnoses, and identity/paternitA testing -- assuming that the 
local surrounding DNA sequence is Woown. 

^ While the invention has been described in connection with 
specific embodiments thereof, it will be understood that it is 
capable of further modifications and this application is intended to 
cover any variations, uses, or adaptations of the invention 
following, in general, the principles of the invention and including 
such departures from the present disclosure as come within known 
or customary practice within the art to which the invention 
pertains and as may be applied to the essential features 
hereinbefore set forth and as follows in the scope of the appended 
claims. 



